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The low-temperature Hall resistivity p^y of La2/3Ai/3Mn03 single crystals (where A stands for 
Ca, Pb and Ca, or Sr) can be separated into Ordinary and Anomalous contributions, giving rise to 
Ordinary and Anomalous Hall effects, respectively. However, no such decomposition is possible near 
the Curie temperature which, in these systems, is close to metal-to-insulator transition. Rather, 
for all of these compounds and to a good approximation, the pxy data at various temperatures and 
magnetic fields collapse (up to an overall scale), on to a single function of the reduced magnetization 
m = M/Msat, the extremum of this function lying at m « 0.4. A new mechanism for the Anomalous 
Hall Effect in the inelastic hopping regime, which reproduces these scaling curves, is identified. This 
mechanism, which is an extension of Holstein's model for the Ordinary Hall effect in the hopping 
regime, arises from the combined effects of the double-exchange-induced quantal phase in triads of 
Mn ions and spin-orbit interactions. We identify processes that lead to the Anomalous Hall Effect 
for localized carriers and, along the way, analyze issues of quantum interference in the presence of 
phonon-assisted hopping. Our results suggest that, near the ferromagnet-to-paramagnet transition, 
it is appropriate to describe transport in manganites in terms of carrier hopping between states 
that are localized due to combined effect of magnetic and non-magnetic disorder. We attribute 
the qualitative variations in resistivity characteristics across manganite compounds to the differing 
strengths of their carrier self-trapping, and conclude that both disorder-induced localization and 
self-trapping effects are important for transport. 

PACS No: 75.30.Vn, 72. 20. My, 7f .38.-fi, 03.65. Bz, 7f .23.An 



I. INTRODUCTION AND OVERVIEW 

Numerous recent studies have focused on the Hall 
effect in the family of doped manganese oxides 
Lai_3;Aj:Mn03 (in which A stands for Ca, Sr or Pb), 
famous for its colossal magnetoresistance (CMR) [np] 
and accompanying ferromagnet-to-paramagne t (FP) and 
metal-insulator (MI) transitions |,|,|,|j7l|,|,|ig,y,|| . 
Despite substantial variations in, e.g., the ferromagnet- 
to-paramagnet transition temperature Tc and residual 
resistivity across this manganite family, measurements 
of the Hall effect reveal unusual features in both their 
metallic and insulating regimes. An example of the Hall 
effect data is shown in Fig. |l|. In the metallic state the 
Hall (i.e. transverse) resistivity pxy at lowest tempera- 
tures (curve at 10 K in Fig. |^) exhibit just the ordinary 
Hall effect (OHE), proportional to the external magnetic 
field B. At higher temperatures in the metallic phase, 
the Hall resistivity can be separated into the sum of: 
(i) a (positive) ordinary Hall effect, and (ii) a (negative) 
anomalous Hall effect (AHE), proportional to the mag- 
netization M, as shown for the curve at 200K on Fig. ||. 
The effective density of carrier holes, as deduced from 
the slope of OH resistivity, is typically found to be sev- 
eral times larger than that set by the nominal doping 
level. 



This difference has been attributed to the effects 
of charge compensation and Fermi-surface shape 
The AHE is commonly observed in ferromagnets, but 
the sign and the magnitude of the AHE in man- 
ganites stand in contradiction to conventional theories 
based on skew-scattering [^3|,^ 15 1^ or side-jump pro- 
cesses |]l5| , |l7| . p"8Hl9| . Most striking is the rapid increase 
in the prominence of the AHE that occurs at temper- 
atures T close to Tc- In this range of temperatures, 
Pxy can no longer be simply separated into ordinary and 
anomalous parts, as can be seen from the curve at 300 K 
in Fig. |l} For temperatures well above Tc, Pxy again 
becomes linear in B, although its sign is now negative 
SliiiSS- The corresponding Hall coefficient Rh 
(= Pxy/ B) decreases exponentially with increasing tem- 
perature in this regime, and a previous study identified 
a clear crossover from non-polaronic to polaronic charge 
transport at around 1.4 Tc [Q. 

The purpose of the present Paper is to address the issue 
of charge-carrier motion in manganites, both experimen- 
tally and theoretically, focusing on the vicinity of the FP 
and MI transitions, from the vantage point afforded by 
the Hall effect. Our experimental results have led us to 
focus on the anomalous contribution to the Hall effect, 
and to develop a microscopic theoretical picture of the 
charge-carrier motion that gives rise to this contribution 
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in manganites. The picture that emerges is one in which 
the essential character of charge-carrier motion is inelas- 
tic hopping between states localized due to magnetic and 
other sources of disorder. 
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FIG. 1. Hall resistivity of manganites versus magnetic field 
for a selection of temperatures. At 10 K the Hall effect is or- 
dinary; the slope extrapolates to the origin. At 200 K the 
Hall effect has both ordinary and anomalous components; 
the slope does not extrapolate to the origin, the offset sig- 
naling the anomalous Hall effect. At 300 K it is not simple 
to separate the Hall resistivity into ordinary and anomalous 
components. 

In order to explain the Hall effect in the manganites in 
the vicinity of Tc, it is necessary to understand how the 
nature of the charge-carrier states are influenced by the 
magnetic order of the system. In this regard, the double- 
exchange interaction (DEI), which makes charge-carrier 
motion of Mn outer-shell carriers sensitive to magnetic 
alignment of core 3/2 spins of Mn ions (Hund rules lead to 
alignment of spins in three inner orbitals resulting in core 
spin 3/2), has long been known to play a key role in trans- 
port in manganites, having been introduced b y Z ener ]2C| ] 
and elaborated by Anderson and Hasegawa tel and De 
Gennes p^ ]. Therefore, our approach to exploration of 
the anomalous Hall effect in manganites is based on the 
picture of hopping conduction in the presence of double- 
exchange interaction. 

A picture of the ordinary Hall effect in hopping con- 
ductors was developed long ago by Holstein ||23| , in which 
the critical ingredient is the Aharonov-Bohm quantal 
phase acquired as charge-carriers hop in the pres- 
ence of magnetic field around closed sequences of local- 
ized states. The theoretical work reported here amounts 
to the generalization of Holstein's ideas suited to DE sys- 
tems. The primary distinctions from Holstein's ideas 
are as follows: (i) Localization is now, to a great ex- 
tent, caused by magnetic disorder in the orientation of 



core spins (and the attendant randomization of hopping 
amplitudes); the effects of magnetic disorder are facili- 
tated by static disorder, and accompanied by polaronic 
effects, (ii) The relevant quantum-mechanical phases 
now arise via the quantal version of the Pancharatnam 
phase [|25| , ^ . (iii) In order for a net Hall effect to result, 
account must be taken of Dzyaloshinskii-Moriya spin- 
orbit coupling |27|. The AHE mechanism that we pro- 
pose arises in hopping regime in systems with localized 
states, and is the only possible microscopic mechanism of 
AHE in such systems. A brief account of this work was 
published in Refs. 0,|ll|- 

From the perspective of symmetry, it is 
well-known that spin-orbit interactions lead to 



AHE Pj,y,|a,h«,h7 18:1J. The appreciation that a 



spin-generated geometric phase, in addition to spin-orbit 
interactions, is an essential ingredient for the devel- 
opment of a theory of the AHE in DE systems dates 
back to the 1998 manuscript of Kim, Majumdar, Millis 
and Shraiman p8[ |. From the perspective of the mi- 
croscopic mechanism of the AHE, both Ref. and 
the paper of Ye et al. that superceded it, invoke 
a field-theoretic scheme for integrating out the charge- 
carrier motion and, therefore, were intended for metallic 
regimes (i.e. regimes in which charge transport occurs via 
delocalized states). By contrast, the present work con- 
siders nonmetallic regimes (i.e. regimes in which charge 
transport occurs via inelastically-assited hopping be- 
tween localized states). Elsewhere, we shall address the 
issue of the microscopic mechanism of the AHE in the 
metallic regime [ pO[ . 

The microscopic mechanism of AHE that we propose 
for systems with localized states leads to remarkable pre- 
diction: the Hall resistivity p^y depends on the tempera- 
ture and the magnetic field solely through the magnetiza- 
tion M, i.e., pxy — pxy{M{H,T)). This universal scahng 
has been observed experimentally in manganites. Here, 
we provide a detailed discussion of our theoretical picture 
of hopping transport in manganites. Further, the present 
paper reports on measurements made on additional com- 
pounds having lower and higher transition temperatures 
and provides an analysis of these data in terms of our 
theoretical picture The universal scaling rela- 

tion between p^y and M reported for is shown to hold 
for the manganese oxides La2/3Cai/3Mn03 (LCMO), 
Laa/sSri/aMnOa (LSMO), and Laa/aCai/gPbz/gMnOa 
(LPMO). Although data on the Hall effect in these com- 
pounds have universal features, the the temperature de- 
pendence of resistivity in LSMO is different from that in 
LCMO and LPMO. This behavior is due to different size 
of dopant ions which results in different static disorder, 
different carrier localization length, and, accordingly, dif- 
ferent strength of self-trapping due to lattice effects. We 
believe that the accuracy of the results concerning the 
AHE based on inelastic charge-carrier hopping between 
states localized due to magnetic and nonmagnetic sources 
of disorder suggests that the dominant mechanism for 
charge transport in the transition regime is indeed inelas- 
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tic charge-carrier hopping between locahzed states, which 
differs quaUtatively from the picture of metaUic conduc- 
tion perturbed by double-exchange interactions. As for 
polaronic effects, depending on the compound, they may 
set in soon as localization length is of the order of lattice 
constant. These effects (or their absence) are crucial for 
the character of the temperature behavior in the range 
of high temperatures above the FP and MI transitions. 
Polaronic effects do not affect the universal scaling of the 
Anomalous Hall resistivity. At the same time, scaling 
of the AH resistivity of the type observed in the CMR 
regime is not contained in conventional models of the 
AHE in the metals (i.e., those based on skew-scattering 
and side-jump mechanisms). Neither is this scaling con- 
tained in a Berry phase mechanism for the AHE in the 
metallic phase, discussed in Refs. p8| , p9[ . We regard this 
as further evidence against the viability of any metallic- 
based picture of transport in the transition regime. 

How does the present work relate to earlier work on 
charge transport in the CMR regime? Attempting to 
build on the early key insight that DE plays a central 
role, Millis et al. considered transport in DE systems 
within the framework of the coherent potential approxi- 
mation (CPA). Making the CPA in the present context 
amounts to replacing the charge and magnetic system by 
an effective one, involving only the charge subsystem, in 
which the conduction band width depends on the magne- 
tization but there is no other effect of the magnetic sector. 
Thus, any resistivity obtained via such a scheme is simply 
whatever the resistivity of the charge sector was, reduced 
by an extent that depends on the magnetic order via a 
renormalization of the bandwidth. The picture enforced 
by this approach is that the fundamental mechanism for 
charge transport is metallic in nature. Naturally, the 
CPA approach is unable to yield a colossal magne- 
toresistance, although it can provide factors on the order 
of unity. What it explicitly omits is any resistivity mech- 
anism arising from localization due to magnetic disorder, 
as noted by Varma |32|. Rather than appeal to such a 
localization process, Millis et al. ||3^ proposed that the 
magnetization-dependent reduction of the bandwidth in- 
vites lattice effects. Specifically, the (now magnetically) 
heavy charge carriers would be more susceptible to self- 
trapping by a large Jahn- Teller lattice distortion, which 
would cause a metal-insulator transition via polaronic 
collapse of the conduction bandwidth. 

Accepting, for the moment, the notion that charge 
transport in the transition regime is indeed accomplished 
by lattice polarons, let us ask what the consequences 
would be for the resistivity. According to theory of 
polaronic transport, developed in series of papers in 
1960's by Holstein alone [Q, with Friedman |35| and 
with Emin | |3^ , polaronic- type conduction manifests it- 
self via a specific temperature-dependence of the longitu- 
dinal and Hall resistivities, being activated in character 
with a definite relationship between the activation con- 
stants for these resistivities. Following the proposal of 
polaronic-type conduction by Millis et al., experimental 



tests of these temperature dependences were performed. 
Initial results Q in LPMO over the range of tempera- 
tures high above the MI and FP transitions seemed in 
accordance with the polaronic picture. However, recent 
extensive measurements at lower temperatures, in transi- 
tion regime fl^ ], demonstrate that, at least in this regime, 
the temperature dependence of the longitudinal and Hall 
resistivities cannot be explained in terms of polaronic 
picture alone. Furthermore, even at high temperatures, 
polaron-based picture is not compatible with experimen- 
tal data for LSMO samples. 

With the pictures of charge transport in the CMR 
regime based on either the magnetization-dependent re- 
duction of the bandwidth or on polarons alone invali- 
dated, what remains is the possibility of constructing 
a valid picture based on non-polaronic localization of 
charge carriers. Strong evidence supporting such a pic- 
ture comes from a simple estimate of the scattering 
time (i.e. the scattering-induced conduction-band broad- 
ening), which indicates that, in the transition regime, 
the band broadening exceeds the band width (i.e. the 
mean free path is shorter than the Fermi wavelength) 
so that the resistivity exceeds the Mott-Ioffe-Regel limit 
and, hence, the conduction cannot be metallic. There- 
fore, one needs to search for insulating transport mech- 
anisms and, specifically, the origins of carrier localiza- 
tion that are distinct from polaronic effects. (We note 
that in compounds in which the Jahn- Teller distortion 
is not symmetry-allowed, this is especially important). 
Such localization can result from both magnetic disorder 
(i.e. due to lack of core-spin alignment) and non-magnetic 
disorder (i.e. static potential disorder due, e.g., to dop- 
ing) H- 

While the localizing influence of non-magnetic disor- 
der on c harg e transport has been thoroughly investi- 
gated p8| , ^9|] , the influence of magnetic disorder i s less 
well known, and we shall discuss it in detail in Sec. IIIB. 



For now, we simply mention that the magnetic disorder 
in core 3/2 spin orientation experienced by the outer-shell 
charge carriers arises via the DE interaction from fluctua- 
tions around the ferromagnetic state that build up as the 
FP transition is approached from the low temperature 
side. Of course, these fluctuations are dynamical, but 
they are slow, compared with characteristic timescales 
for outer-shell charge-carrier motion. Thus, for the pur- 
poses of analyzing the influence of the magnetic sector on 
charge transport, it is appropriate to regard orientations 
of the core spins on the Mn ions as quenched variables. 
The resulting magnetic disorder takes the form of ran- 
domness in the off-diagonal hopping matrix elements for 
the charge carriers. By contrast, nonmagnetic disorder 
occurs due to randomness in the substitution of La by 
dopant ions (e.g. Sr, Pb or Ca), and gives rise to the 
more familiar diagonal (Anderson-type) disorder. Elec- 
tronic states in systems with off-diagonal disorder were 
first considered by Lifshitz |Q, who showed that local- 
ized states arise in the band tail. The physical picture of 
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carrier states in manganites must encompass both mag- 
netic and nonmagnetic disorder, possibly facilitated by 
Coulomb effects, which (jointly or severally) can result 
in carrier localization. 

If carriers are localized then they can still participate 
in transport, but it is by hopping from one localized state 
to another, assisted by one or more inelastic agents (such 
as phonons). In this case, the longitudinal resistivity is 
determined by the rate of inelastic hopping between oc- 
cupied and unoccupied states |^,^. When carrier lo- 
calization has occurred, and the localization length is of 
order of lattice constant, electronic interaction with lat- 
tice and self-trapping effects can become essential, so that 
at high temperature resistivity is determined by small po- 
larons. However, transitional regime is greatly affected 
by carrier localization of non-polaronic origin. We no- 
tice that in general one should distinguish between Jahn- 
Teller polarons and Holstein breathing mode polarons: 
the presence of the former depends on symmetry of the 
system, the latter arise independent of the underlying 
symmetry. In manganites, both types of polarons are 
capable of facilitating carrier localization by magnetic 
and non-magnetic disorder; when carriers are localized 
on lattice constant scale, Holstein polarons govern the 
temperature dependence of resistivity deep in the insu- 
lating phase. 

In the present Paper, we shall not consider metallic 
manganites, and restrict our consideration to the inelas- 
tic hopping transport regime. The discussion of the Hall 
effect in metallic ferromagnets will be presented else- 
where . The present Paper is organized as follow s. In 
we describe the experimental setup and in Sec. [II Ij| 
we present experimental data on the longitudinal resistiv- 
ity, magnetization and Hall resistivity in three different 
manganite compounds having distinct transition temper- 
atures. Section [II is organized into several subsections, 
in which we describe models of disorder in manganites, 
issues related to the localization of carriers and the hop- 
ping transport mechanism, as well as the quantal Pan- 
charatnam phase, spin-orbit Dzyaloshinski-Moriya inter- 
actions, and the universal scaling of the Hall resistivity. 
In Sec. IV we discuss the correspondence between our 



theoretical and experimental results. 



II. EXPERIMENTS ON TRANSPORT AND 
MAGNETIC PROPERTIES OF MANGANITES 

A. Experimental method 

In the experimental part of this study, single crys- 
tals of various manganites were used in order to avoid 
extrinsic effects from grain boundaries or strains. In 
single crystals, simultaneous measurements of trans- 
port and magnetic properties permits us to find a pre- 
cise dependence of transport coefficients on the sam- 
ple magnetization. We have measured the longitu- 



dinal and Hall resistivities and the magnetization of 
three different crystals with various transition temper- 
atures. Lao.rCao.sMnOa (LCMO) and Lao.TSro.aMnOa 
(LSMO) single crystals were prepared by the floating- 
zone method. Lao.67(Ca,Pb)o.33Mn03 (LPMO) single 
crystals were grown from 50/50 PbF2/PbO flux. More 
details on the sampl e grow th and basic properties can 
be found elsewhere |42|j43|| . All specimens used in the 
measurements were cut along crystalline axes into bar 
shapes from larger pre-oriented crystals. Contact pads 
for Hall resistivity measurements were made by sputter- 
ing «1500 A of gold through a mask. Gold wires of 50 
/im diameter are then attached using slowly drying silver 
paints. Typical contact resistances after annealing were 
about 1 Q ai room temperature. We adopted a low- 
frequency (39 Hz) ac method for the measurements. The 
transverse voltage signal was flrst nulled at zero field at 
each temperature below 400 K by a potentiometer, and 
the change in the transverse voltage was recorded as H 
was swept from 4-7 T to -7 T and back for averaging. 
Following the transport measurements, sample magneti- 
zations were measured by a 7 T SQUID magnetometer 
on the same samples. 



B. Experimental Results 



Figure]^ shows the temperature dependences of magne- 
tization measured at 1 T and 7 T. All three samples show 
ferromagnetic-to-paramagnetic phase transitions. The 
Curie temperatures Tc were determined by scaling anal- 
ysis on high field M{H) curves near the transition, and 
the results are shown in Table I. As Tc decreases, the 
transition becomes sharper, resulting in anomalous crit- 
ical exponents |Q. 

The temperature dependences of the longitudinal resis- 
tivities pxx for the same set of samples under zero mag- 
netic field and under 7 T are shown in Fig. ^. LCMO 
and LPMO show metal- insulator transitions near Tp, 
whereas LSMO shows an inflection at T^, but p^x con- 
tinues to increase with increasing temperature above Tq- 
The metal-insulator transition temperatures Tmi , deter- 
mined by the maximum in the rate of change in the 
temperature dependence of the longitudinal resistivity 
dpxx/dT under zero magnetic field, were slightly higher 
than corresponding Tc's (Table I). Also shown in Ta- 
ble I are pxx minima (occurring at the lowest tempera- 
tures) and magnetoresistivity (MR) maxima [defined by 
{pxxiO T)-Pxxi7 T))/pxxi7 T)]. The observed decrease 
in Tc correlates with the overall increases of resistivity 
and MR, which can be clearly seen in Fig. ||. 

Despite differences in Tc, Tmi and temperature de- 
pendence of longitudinal resistivity across the three com- 
pounds, the Hall resistivity of compounds with doping 
that corresponds to maximal Tc show similar tempera- 
ture and field dependences, as shown in Figs. ^, ^ and||. 
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At low temperatures, pxy is positive and linear in mag- and negligible anomalous Hall contribution, 
netic field, the sign indicating hole-like charge carriers, 
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Lao.rCao.sMnOa 

Lao.67(Ca,Pb)o.33Mn03 

Lao.rSro.aMnOa 


216.2 K 
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287.5 K 
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FIG. 2. Temperature dependences of normalized magneti- 
zation M/MsiA under 1 T (solid lines) and under 7 T (dotted 
lines) magnetic fields. 
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FIG. 3. Temperature dependences of longitudinal resistiv- 
ity Pxi, under zero magnetic field (solid lines) and under 7T 
(dotted lines). 

As the temperature is increased, the high-field slope is 
roughly the same. However, the increasing, negative, 
contribution to Pxy shifts it downward. Quantitatively, 
Pxy can be expressed as a sum of an ordinary contribution 
parametrized by i?o (T') and an anomalous contribution 
parametrized by i?s (r) |4|] : 



Pxy {B, T) = i?o (T) B + ^loRs (T) M (H, T) , (2.1) 

where B = /xg [i? + (1 — N)M] , H is the external magnetic 
field, and iV « 1 is the demagnetization factor. As has 
commonly been observed in manganite crystals and thin 
films, the effective charge-carrier density n^ff = I/bRq 
is scattered between 1.0 and 2.4 holes/Mn (see Table I), 
which is much larger than the nominal doping level (of 
0.3-0.33 holes/Mn), presumably due to the effects of the 
anisotropy of the Fermi surface |^]. We refer to our pre- 
vious publications for the discussion on the low temper- 
ature OHE g. 

On further increase of the temperature through Tc, 
Pxy becomes much larger, strongly curving with mag- 
netic field, and the positive, high-magnetic-field contri- 
bution, linear in the field, which would arise from the or- 
dinary Hall effect in a metallic phase due to the Lorentz 
force acting on charge-carriers, disappears. Owing to its 
low Tc, for the LCMO sample we were able to explore 
temperatures far above Tc, where pxy shows a negative 
Hall coefficient, despite the doping of the material being 
by holes. In this range of temperatures the Hall coeffi- 
cient Rh = Pxy{B)/B exhibits activated behavior, with a 
characteristic energy Eh ~ where E„ is the activa- 
tion energy for ordinary conductivity Gxx ■ Similar experi- 
mental results have been obtained in Refs. . In these 
works, investigations of manganite carrier-transport deep 
in the insulating phase have shown that the sign and tem- 
perature dependence of the high-temperature (i.e. above 
1.4rc') Hall coefficient Rh (= Pxy/B) can be explained 
in terms of the adiabatic hopping of small polarons. 

Initially Q, high-temperature transport picture due 
polaron hoppingto was believed to support the proposal 
by Millis et al. ||31|] that the Jahn- Teller distortion which, 
according to symmetry considerations can occur in the 
MnOg octahedra of LaMnOa, is responsible for the insu- 
lating behavior of doped Lai-jjA^jMnOa systems. Hall 
resistivity measurements should be capable of provid- 
ing key evidence for or against the polaronic picture of 
charge transport. According to the theory of this pic- 
ture, the adiabatic hopping of small polarons |^ leads 
to an activation energy Eh characterizing Rh that is 2/3 
of the activation energy Ep characterizing pxx | |3^ , as is 
observed at high temperatures However, recent 
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Hall resistivity measurements, extending into the transi- 
tion region show that the activation energy changes 
abruptly at a cross-over temperature lATc, from the 
polaronic value of '^Ea to a much larger value, 1.7 E^r- 
This clearly marks the breakdown of the small polaron 
picture of charge transport, as shown in Fig. ^ (inset). 
In fact, the effective activation energy of the conduc- 
tivity begins to decrease from a value of Ea- at roughly 
the same cross-over temperature, making even greater 
the discrepancy between the experimental data and the 
small-polaron hopping picture. Even more dramatically, 
the product of the Hall mobility fj,H and the temperature, 
viz., (XhT — —axxRiiT which, according to the small 
polaron picture, should decrease monotonically with de- 
creasing temperature, in fact is found to exhibit a mini- 
mum at the same cross-over temperature. (We shall later 
show that, near Tc, Pxy is determined solely by the sam- 
ple magnetization in all three compounds.) Thus, exper- 
iments in transition region lead to the conclusion that, 
while small polarons are an essential part of the physics of 
transport in manganites at high temperatures, they can- 
not provide a complete picture of the metal-to-insulator 
transition. 
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FIG. 4. Main panel: Hall resistivity p^y of LCMO as a 
function of magnetic field at the indicated temperatures. In- 
set: Activated behavior of the high temperature Hall coeffi- 
cient Rh- 

More generally, as discussed by Varma |Q, there ex- 
ist double-exchange systems, such as TmSea;Tei_3;, in 
which transport phenomena observed in manganese ox- 
ides are also observed but Jahn- Teller distortions, lead- 
ing to small polarons, are not symmetry-allowed. At the 
same time, if the carrier localization length becomes of 
order of lattice constant, lattice effects in the form of Hol- 
stein breathing mode polarons arise naturally. This al- 
lows to explain why the high-temperature regime in some 
of manganese compounds exhibits longitudinal and Hall 
resistivities characterized by thermally-activated behav- 
ior which is qualitatively and quantitatively consistent 
with that caused by polaronic transport mechanism. 



0.0 



LPMO 



-0.5 - 




33D K 



1234567B 
S(T) 

FIG. 5. Hall resistivity p^y of LPMO as a function of mag- 
netic field at the indicated temperatures. 
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Before turning to the theoretical picture of transport 
in manganites and, specifically, of AHE for localized 
carriers, we pause to examine whether the theoretical 
model of the AHE proposed by Ye et al. which is 
based on a metallic view of charge transport, is consis- 
tent with our experimental data. In that model, follow- 
ing an earlier model by Kim et al it is assumed 
that tight-binding charge-carriers propagate coherently 
through a smoothly varying magnetization texture which 
has the effect of introducing a Berry phase gauge poten- 
tial |4^,^,||,||,|o|. A central prediction of the model 
due to Ye et al. is that a peak should occur in Rs{T) 
above Tc, along with a singularity in the slope at Tc, 
i.e., dRs/dT ~ |1 - T/Td"" + c, where a is the spe- 
cific heat exponent. To test this prediction, we measured 
the low magnetic field (< 0.5 T) magnetization and Hall 
resistivity of our most metallic sample, LSMO near Tc 
(see Fig. From the behavior of p^y and M in the zero- 
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field limit, we determined Rs = [dp^v/dM). In constrast 
to the prior report by Matl et al. Q, in this "metal-to- 
metal" transition system we do not find Rs to be pro- 
portional to pxx- As seen on Fig. p^x flattens at the 
temperature at which the resistive transition is complete 
(i.e. — 368 rniK). This allows us to conclude that 
neither a constant Rs/ Pxx nor a peak in Rs is a common 
feature in manganites. We note that close to T^, it is pos- 
sible to express Rs(T) as a power law, {l-T/T^)°-^^ + A 
(see Fig. 0, inset). However above Tq = 368K, which is 
significantly higher than both Tc and Tmi (see Table I) , 
Rs is constant, and the fit does not correspond to the 
inflection point predicted in Ref. | |29| ]. 
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FIG. 7. Main panel: The anomalous Hall coefficient Rs 
(symbols) compared with the longitudinal resistivity pxx 
(solid line). Inset: The critical behavior of Rs. 



III. THEORY OF HOPPING 
MAGNETOTRANSPORT IN MANGANITES 

The aim of the present section is to develop a pic- 
ture of the Hall effect in manganites, to test this picture 
through comparison with experimental data and, hence, 
to build as completely as possible a general picture of the 
of charge transport in manganites in the ferromagnet- 
to-paramagnet transition regime. Among the results we 
shall obtain, perhaps the most striking is the universal 
scaling of the magnetization-dependent Hall resistivity, 
which we explain should hold in the regime where charge 
transport proceeds primarily via hopping between local- 
ized states. Such universal scaling has been observed 
experimentally. 



A. Disorder and interactions in manganites 

Manganites are extremely complicated materials, and 
a bewildering variety of behaviors occurs in them, as 
the doping level, temperature or magnetic field is var- 
ied. Here, we focus on those manganites that exhibit a 
transition from a ferromagnetic metal to a paramagnetic 



insulator, controlled by temperature, as occurs in man- 
ganite compounds, doped with Ca or Sr or (Pb and Ca) 
substituting for La, at doping levels of around 1/3. This 
doping, of course, results in several sources of static disor- 
der: (i) the dopant ions subsitute randomly for La; and 
(ii) the lattice distortion around the two ionizations of 
Mn (viz. Mn'^+ and Mn''+) is distinct (i.e. the breathing- 
mode effect). This disorder leads to local variations in the 
amplitudes for the hopping processes that carry charges 
between magnetic ions. Furthermore, (iii) any clustering 
of dopants in the randomly doped lattice would lead to 
fluctuations in the carrier density. These sources produce 
nonmagnetic disorder. 

Along with the motion of charge carriers there is also 
the motion of core spins. In the present context, we 
believe that it is profitable to treat theses core spins 
classically, and to regard the carrier dynamics as being 
much faster than the spin dynamics, so that the carrier 
motion can be pictured as taking place within a frozen 
core-spin configuration that is randomized owing to ther- 
mal fiuctuations (i.e. we adopt a quasi-static approach). 
We term such random magnetic configurations "magnetic 
disorder'.' We note that spin-spin correlation times have 
been obtained experimantally from muon spin relaxation 
and neutron spin echo data Q ; these experiments show 
that, indeed, spin dynamics is slow. 

Strong thermal fiuctuations render typical instanta- 
neous configurations of the spins rather inhomogeneous. 
Among these fluctuations, there are the "hedgehog" 
excitations which, owing to their topological stabil- 
ity, are long-lived, and become more numerous as the 
ferromagnet-to-paramagnet transition is approached . 
Due to the resulting magnetic inhomogeneity, the carrier- 
hopping matrix elements are reduced. 

As we shall discuss in the following subsection, the 
presence of nonmagnetic and magnetic disorder both sup- 
port the notion that the carrier states are localized at 
temperatures near to the (zero-magnetic-field) FP tran- 
sition, as well as at higher temperatures. Such localiza- 
tion of carriers can explain the resistive transition which, 
in turn, leads to the disappearance of double-exchange- 
induced ferromagnetic spin-correlations, at least on spa- 
tial scales larger than the localization length. We note 
that in a series of papers |53| Furukawa has considered 
the issues of carrier states and transport in manganites 
by using a dynamical generalization of the coherent po- 
tential approximation, arriving at the conclusion that the 
transport properties of manganites can be explained in 
terms of the scattering of metallic carriers by magnetic 
randomness. We, however, believe that the fact that (in 
the range of temperatures marking the transition regime) 
the resistivity exceeds the Mott-Ioffe-Regel limit renders 
any approach founded on the scattering of delocalized 
carrier-states to be inconsistent. 

As discussed above in Sec. |, localization effects related 
to Jahn- Teller distortions and small-polaron formation 
cannot explain certain central experimental data in the 
transition regime, including the temperature-dependence 
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of both the resistivity and the Hall effect. Therefore, one 
is forced to consider alternative mechanisms that can lead 
to the breakdown of metallic conductivity and can also 
serve as an origin of various universal transport proper- 
ties that have been observed in double-exchange systems, 
including, e.g., those in which Jahn- Teller distortions are 
symmetry- forbidden. For these reasons, we now give a 
discussion of the physics of disorder-induced carrier lo- 
calization in manganites. 



B. Disorder-induced carrier localization in 
manganites 

To see why it is useful to regard charge transport as 
taking place in a frozen random background of core-spin 
orientations, let us imagine the spin configuration to be 
truly static. In the transition regime, a typical spin con- 
figuration is rather inhomogeneous and, hence, we expect 
carriers to be localized. Support for this notion comes 
from the close similarity between transport in mangan- 
ites and systems of randomly located identical impuri- 
ties (i.e. off-diagonal disorder). For the latter, localiza- 
tion has been established via the work of Lifshitz pO[ . 
Although spin-induced randomness in manganites [aris- 
ing from the random double-exchange factors of cos(6'/2), 
where i s the angle between core spins on Mn ions, see, 
e.g.. Sec . |III E| is weaker than the randomness considered 
in Ref. ||4Cl|], we expect the two systems to exhibit qual- 
itatively similar localization behavior. Furthermore, the 
condition for localization in the Lifshitz model (viz. that 
the characteristic spatial scale of the outer-shell wave- 
functions in isolated Mn ions be much smaller than dis- 
tance between sites) is well obeyed in manganites. There- 
fore, provided that there is appreciable randomness in 
the core-spin orientations, transport properties should 
be determined by the short-distance physics of clusters 
of ions and by magnetic correlations between such clus- 
ters. Moreover, nonmagnetic disorder and possible states 
bound to the subsituting A-site ions are capable of am- 
plifying the trend towards localization [ ^2[|37| ]. In Fig. ^ 
we present a one-dimentional caricature of disorder in 
manganites, in which diagonal and off-diagonal disor- 
der coexist. Sheng et al. included both magnetic and 
non-magnetic disorder and applied one-parameter scal- 
ing theory |54j and finite-size scaling ideas |^ in order 
to investigate carrier localization in manganites numeri- 
cally. Sheng et al. found that, in the presence of magnetic 
disorder, an Anderson metal-insulator transition accom- 
panies the ferromagnet-to-paramagnet transition. They 
also observed an interesting correlation between Tc and 
the residual resistivity, which is determined by nonmag- 
netic disorder, viz., the larger the residual resistivity, the 
smaller the Tc ; this agrees well with the original double- 
exchange picture, in which carrier motion promotes fer- 
romagnetism whereas disorder resists electronic motion 
and, therefore, does not promote ferromagnetism. 



~i_n_r 

FIG. 8. Schematic one-dimensional picture of disorder 
in manganites. (a) Anderson model with on-site disorder; 
(b) Lifshitz model with random ionic locations (leading to 
random transfer integrals); (c) model of disorder in mangan- 
ites, showing on-site disorder and magnetic disorder (leading 
to random transfer integrals) 

To what extent can one regard charge transport as tak- 
ing place in a frozen spin configuration? Equivalently, 
what are the characteristic time scales for magnetic and 
charge dynamics? For charge dynamics the time scale 
is ti/t, where t is a characteristic magnitude of the hop- 
ping matrix element, i.e., the shortest relevant time scale. 
For magnetic dynamics, the issue is more complicated. 
However, the shortest time scale is presumably Ti/fceT, 
which, in the regime of interest to us, is greater than fi/t. 
Furthermore, as discussed by Lau and Dasgupta in 
three-dimensional magnetic systems the ferromagnet-to- 
paramagnet transition involves very long-lived topologi- 
cal excitations (i.e. magnetic hedgehogs) which, as 
we shall see, are particularly significant for the AHE. 
Hence, we see that, to a first approximation, one can re- 
gard charge transport as taking place in a frozen spin 
configuration. Corrections to charge dynamics, due to 
magnetic dynamics (as well as feedback into the mag- 
netic dynamics sector) can, if necessary, be treated by 
going beyond the Born-Oppenheimer approximation. 



C. Percolation- hopping scenario of transport 
phenomena in manganites 

As discussed in the previous subsection, carrier states 
in manganites are effectively localized throughout the 
transition regime. (By effectively localized we mean lo- 
calized on timescales short compared with that required 
for the reconfiguration of the magnetic degrees of free- 
dom.) Therefore, in this regime transport of carriers oc- 
curs via inelastic hopping (i.e. hopping that is assisted by 
some inelastic agent such as a phonon). In the following 
Sec. HI CI, we describe a picture of hopping transport 
app lied to the setting of manganites. Following this, in 
Sec HI C 2 , we use this picture to discuss a scenario for 



transport phenomena in LCMO and LPMO, which are 
materials in which polar onic e ffects are believed to be 
important. Then, in Sec. [II C 3| , we describe the scenario 
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of hopping transport suitable for application to LSMO, a 
material in which signatures of polaronic effects are ab- 
sent. By contrasting LCMO and LPMO with LSMO we 
draw some general conclusions concerning the role played 
by polarons in various manganite compounds. 



1. Hopping transport picture in manganites 



Hopping transport models based on percolation theory 
have been successfully applied to transport in systems 
with states localized by static disorder |3|^^^. A pecu- 
liarity of manganite systems is that the wavefunctions of 
states localized in the vicinity of Mn ions depend on the 
orientations of the Mn core spins on these ions. Inelas- 
tic agents (e.g. phonons) lead to hopping between these 
localized states, the amplitude for such hopping being 
characterized by matrix elements of the carricr-phonon 
interaction. Therefore, hopping probabilities and rates 
are determined by the orientations of the core spins. In 
the d ouble -exchange picture [discussed in more detail in 
Sec. ([HE)], the rate Wij of carrier- hopping between ions 
i and j, whose core spins (which we treat classically) form 
an angle 6, is proportional to cos^(0/2). 
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FIG. 9. Schematic picture of the conducting network. 
Zig-zag line connecting sites i and j represents the resistivity 
of the bond between these sites. 

As in the standard percolative approach to transport 
in impurity systems (i.e. the Miller-Abrahams resistive- 
network approach |^,^,^), one connects every nearby 
pair of Mn ions by a bond, and assigns a resistivity to 
each of these bonds. The charge current Jjj between ions 
i and j is then given by: 



to linear order in the electric field. Hence, one can ob- 
tain (see e.g., [^) an expression for the electric current 
in Ohmic form: 



kT 

- 



(3.2a) 
(3.2b) 



where Ui — Uj is the potential difference between sites i 
and j in the presence of the electric field, and is the 
(E — 0) transition rate. Therefore, the resistance Rij 
of the bond between ions i and j (cf. Fig. H) is propor- 
tional to 1/Wij, where we have, for the sake of brevity, 
omitted the superscript (which indicates zero electric 
fields quantities) . The resistances of the bonds constitute 
a resistive network on which carriers move by taking the 
least resistive paths, i.e., the transport is percolative in 
character. The conductivity of the sample is entirely de- 
termined by a set of hopping resistivities Rij . 



2. Scenario of transport properties in manganites 

In the hopping regime, the following scenario for trans- 
port properties of manganites can be envisioned: 




J^i = e (W,j 



(3.1) 



7(K) 

FIG. 10. Transport regimes in manganites exeplified by 
LCMO. I: High-temperature regime, in which transport is via 
polarons. II: Crossover to inelastic hopping of charge-carriers 
between localized states. Ill: Maximal resitivity. In this 
regime the loss of inelastic agents is compensated by the 
growth of magnetic order and, hence, the emergence of a con- 
ducting network (see text for details). IV: Rapid growth of 
the conducting network. V: Saturation of the conducting net- 
work. VI: Crossover to the metallic regime. 



In the presence of an applied electric field E and for a 
closed external circuit the system out of equilibrium and 
the charge current is nonzero. If the electric field is suffi- 
ciently small (i.e. eE • (Rj — Ri) <C kT), one can expand 
the charge carrier energies and site occupation numbers 



(i) In the paramagnetic insulating state, i.e., in region I 
on Fig. 1^, the percolative inelastic motion of strongly 
localized carriers is suppressed by magnetic randomness, 
via the DE interaction. Although there may exist small 
clusters of magnetic ions with spins aligned in some 
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direction, neighboring spins (or the spins of neighbor- 
ing aligned smaU clusters) are weakly correlated, and 
are therefore predominantly have a large angle between 
them. Thus, the resistance of the corresponding resistive 
bonds is generally large. Hence, the clusters are isolated 
from each other (i.e. no percolative path exists for which 
core spins on neighboring Mn ions are approximately 
aligned), so that outer-shell carriers cannot hop along any 
path of bonds without encountering a large resistance. 
Furthermore, if the localization length is on the order of 
the lattice constant then carrier interactions with the lat- 
tice and carrier self-trapping via small-polaron formation 
become important. (Below, we shall discuss the situa- 
tion in which the localization length is larger than the 
lattice constant.) Carrier self-trapping, when it occurs, 
does so on very weak (i.e. very resistive) bonds in the re- 
sistive network. Deep in the insulating regime, all paths 
of bonds encounter regions in which carriers are self- 
trapped. Under these conditions, transport occurs via 
the rather infrequent hopping of small polarons, which 
leads to the thermally-activated temperature-dependence 
of the longitudinal resistivity and Hall coefficient. In the 
LCMO and LPCMO compounds, the role of polarons 
in transport accounts for the magnitude of the activa- 
tion energy associated with the resistivity, which is above 
100 meV, and is significantly larger than (any fraction of) 
the hopping amplitude t. 

(ii) With decreasing temperature, the inelastic hopping 
becomes less frequent, so that the resistivity grows. How- 
ever, at these temperatures, percolative paths appear 
that do not encounter regions with self-trapped carri- 
ers. Also, the core spins become more aligned with one 
another, and fairly large clusters exist (that do not fea- 
ture large resistances connecting pairs of Mn ions having 
anti- aligned spins). This regime occurs in the range II of 
temperatures in Fig. 

(iii) With continued decrease in temperature, the resis- 
tivity reaches a maximum (region III in Fig. ^) when 
the core spin orientations become sufficiently correlated 
that a tenuous but infinite conducting network emerges. 
Due to the still-strong magnetic disorder, as well as any 
nonmagnetic disorder, the carrier states are still local- 
ized (and lie in the band tail). The localization length 
is on the order of one to two Mn sublattice units. (It is 
important to note that clusters of spins of size two sublat- 
tice units contain some 20 to 30 spins. The alignment of 
these spins can be obtained by applying magnetic fields 
of order only a few Tesla, leading to colossal magnetore- 
sistance.) 

(iv) With yet further reduction in temperature, the re- 
sistivity decreases abruptly, in line with the standard 
percolation picture [^41|, as more and more inelastic 
hopping-paths become available to carriers, owing to the 
increased alignment of core spins (region IV in Fig. 10). 
The small polarons disappear and, in addition, some 
of what used to be localized states become delocalized, 
so that some carriers now populate the states lying on 
the mobile (i.e. metallic conduction) side of the mobility 



threshold, (v) The abrupt decrease in resistivity with 
temperature slows down as soon as any newly-available 
hopping paths are effectively shunted by the existing con- 
ducting network (region V in Fig. ^|) . 
(vi) Further decrease in temperature leads to further 
core-spin alignment and, ultimately, to a significant den- 
sity of carriers populating states in the conducting part 
of the band and, hence, to the occurrence of the metallic 
state (region VI in Fig. nof). 



3. Scenario of transport properties of manganites: Absence 
and Presence of polarons. 



As described in Sec. |IB|, the temperature dependence 
of the resistivity does not have a universal form across 
all manganite compounds. In particular, in LCMO and 
LPMO, the high-temperature Pxx is thermally activated, 
but in LSMO it is not. It is our opinion that the sce- 
nario described in points (i-vi) in the previous paragraph, 
which involves self-trapping effects due to polaron forma- 
tion, takes place in LCMO and LPMO but not in LSMO. 
This opinion is supported not only by our transport data 
but also by direct neutron-scattering evidence for the co- 
existence of distorted and undistorted Mn-0 octahedra 
in the vicinity of Tc ||5^,^,^ , as well as by the occur- 
rence of a substantial isotope effect in the resistivity, Tc 
and thermal expansion . 

By contrast, LSMO shows no evidence of polarons in 
this suite of experiments. Therefore we propose that in 
LSMO it is only magnetic and static disorder that drive 
the transition between low- and high-resistance states. 
According to the scenario described in the previous para- 
graph, any tendency for the formation of polarons is sup- 
pressed when the localization length at the resistive tran- 
sition (i.e. the inflection point in LSMO) turns out to 
be larger than a few lattice constants. Because of this, 
with further increase in the temperature the localization 
length still has room to decrease due to the suppression 
of ferromagnetic correlations. This would lead to an in- 
crease of the resistivity at temperatures in the immediate 
vicinity above the resistive transition. Such a resistivity 
increase has been observed in LSMO. 

The significance of polarons in LCMO and LPMO, and 
their apparent insignificance in LSMO, is consistent with 
the tendency for polaronic self-trapping to be enhanced 
for reduced A-site ionic radius, as is the case in the 
sequence Sr— >Pb— s-Ca encountered in our experimental 
data. The significance of polarons in LCMO and LPMO 
can also explain why the magnetoresistance of these com- 
pounds is stronger than that of LSMO: In LCMO and 
LPMO, the application of a magnetic field not only re- 
sults in the tendency to delocalize charge carriers by re- 
ducing the magnetic disorder, as it does in LSMO, but 
also such magnetic fields destabilize the polaronic re- 
gions, leading to a more abrupt reconnection of the net- 
work. Furthermore, the presence of polaronic and non- 



10 



polaronic spatial regions in LCMO and LPMO enhances 
disconnection between different parts of the resistive net- 
work, and, due to the DE origin of ferromagnetism in 
these systems, reduces the transition temperature. This 
"boot-strap" collapsing of magnetic order explains why 
the bare double-exchange energy that determined the 
spin-wave dispersion at low temperatures is the same 

even though tLCMO ^ 220K, 



for all manganites [ 

rpLPMO ^ 285K, and T^^mo ^ gggx. We propose that 
the different sizes of the A-site ions which, in particular, 
leads to differing-strengths in the static disorder and dif- 
fering tendencies for self-trapping, is responsible for this 
trend in the transition temperatures. We also note, that 
the existence of polaronic and non-polaronic spatial re- 
gions in LCMO and LPMO explains the success of the 
effective medium approaches in predicting the thermo- 
electric power from the resistivity ||6^ , |6^ , the observa- 
tion of both diffusive and continuum electronic signals 
in Raman scattering [|65| , and the presence of significant 
telegraph noise in the resistivity in these compounds . 
These features are not characteristic for LSMO. 

In Sees. [II F - III I and IV, we shall look at the sce- 



nario that we have just outlined from the vantage point 
afforded by the Hall effect. As, in our opinion, the domi- 
nant transport mechanism in the transition regime occurs 
via inelastic hopping between localized states, regardless 
of whether self-trapping v ia the formation of small po- 
larons occurs, in sections 
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we shall first review the 
basic physical picture underlying the ordinary Hall effect 
in the inelastic hopping regime in systems having poten- 
tial (but not magnetic) disorder. The main issues of this 
discussion are interference of hopping amplitudes in the 
inelastic hopping regime and elucidation of contribution 
to the Hall effect in this regime by using properties of the 
hopping probability with respect to time-reversal sym- 
metry Sees. HID 2 - HID 4. (We obtain expressions for 



hopping amplitudes which differ from those of Holstein 
in ||2^, but this difference is not significant). Readers fa- 
miliar with these issues and Holstein theory of the Hall ef- 
fect can proceed to subsections following this discussion, 
in which we shall provide an extended discussion of our 
picture of the microscopic mechanism of the anomalous 
Hall effect in the hopping regime in manganites, which 
we have recently proposed |10 1^. 



D. Holstein theory of the Hall effect in the hopping 
regime 



of the remaining sites j (which are initially assumed to 
be unoccupied), Wij contains a contribution SWfj that 
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FIG. 11. A triad of sites is the minimal element of the 
conducting network that leads to a Hall electromotive force. 
Owing to the flux $ the carrier current Jh flows perpendic- 
ular to the electric fleld E. The shaded site is occupied; the 
unshaded sites are unoccupied. 

is linear in the applied magnetic field. This dependence 
arises from interference between the amplitude for a di- 
rect (inelastic) transition between the initial and final 
sites and the amplitude for an indirect (inelastic) transi- 
tion, involving the intermediate occupation of the third 
site k. Furthermore, when one applies a magnetic field 
creating a nonzero magnetic flux through the triangle, 
thus introducing an Aharonov-Bohm phase for paths that 
wind around the triangle, the necessity that Wij and Wji 
be equal is lost, even if E _L (R^ — Rj). The Hall current, 
flowing through the bond between sites i and j in a triad 
of sites, in the presence of such electric field, and in mag- 
netic field perpendicular to the plane containing sites i, 
j and fc, as sketched in Fig. |ll|, is given by [c.f. Eq. ([3.1])] 



(3.3) 



Such current would cause an increasing imbalance of pop- 
ulations of sites i and j. However, the balance is restored 
because charge imbalance establishes a chemical poten- 
tial difference Afiij between sites j and i which manifests 



itself as a Hall voltage. Below in Sees. HI F - HI I we shall 



generalize this idea to charge carrier motion in the pres- 
ence of core spins having inhomogeneous orientations. 



1. Ordinary Hall effect: Direct and indirect hopping 



Nearly forty years ago, Holstein |23| observed that to 
capture the ordinary Hall effect in hopping conductors 
requires the analysis of at least triads of sites (i.e. atoms, 
ions, impurities, etc.), and of the attendant Aharonov- 
Bohm (AB) magnetic fiuxes through the polygons whose 
vertices are these sites. What Holstein showed was that 
the probability of the hopping of a charge carrier that is 
initially located on one of three sites i to one or the other 



Let us now address the issue of the Hall effect at an 
elementary level. To do this we consider a system of lo- 
calized carriers, their wavefunctions (r|4'j) (j = 1,2,...) 
being the exact wave functions of the discrete spectrum 
of the electronic Hamiltonian Hd in the presence of ionic 
potentials, potential disorder and magnetic field. (In or- 
der to be specific, we assume, in the present subsection, 
that the carriers are electrons.) 
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Now consider the rates of hopping between these exact 
electronic states caused by the electron-phonon interac- 
tion Wj^k- The necessity of the electron-phonon inter- 
action (or interactions with some other inelastic agent) 
for inducing electron hopping will be discussed below in 
the present subsection. Within the context of hopping 
rates it is valuable to introduce the notion of direct and 
indirect hopping rates. The direct hopping rate Wj!!^;^ is, 
to leading order in the electron-phonon interaction, de- 
termined by the single direct transition amplitude, i.e., 
the electron-phonon interaction matrix element 

C/,, = (vI/,|i7ei-ph|^'fc), (3.4) 

where i?ei-ph is the Hamiltonian of electron-phonon in- 
teraction. In Fermi's Golden Rule approximation W^^i^ 
reads 

WfZk = Y^U,k\'S{E, -E,~ huj), (3.5) 

where {Ei} are the exact energy eigenvalues of Hd, and 
w is a phonon frequency. There are additive corrections 
to the direct transition amplitude, which are associated 
with processes involving phonon-induced scatterings that 
do not change the electronic state. In addition to the di- 
rect transition amplitude, there are amplitudes for indi- 
rect transitions from site j to site fc, which are defined to 
be those amplitudes that involve at least one intermedi- 
ate eigenstate j^';) (but now i is restricted to be neither 
j nor k). Among these, there is a subset of amplitudes 
involving exactly one intermediate eigenstate. Such am- 
plitudes have the following characteristic property: The 
indirect amplitude that proceeds via a third site i nec- 
essarily involves two electron-phonon interaction matrix 
elements Uij and Uki- Direct and indirect transition am- 
plitudes can interfere, and, as we shall describe below, 
lead to the Hall effect. 



Why is it that we need to consider inelastic processes 
when considering the Hall effect? As for the longitudi- 
nal hopping conductivity of localized carriers, it is due to 
electronic quantum transitions between localized carrier 
eigenstates, assisted by phonons (or some other inelas- 
tic agent). The participation of some inelastic agent in 
hopping conduction is required for the following reasons. 
First, owing to carrier localization, the conductivity of 
the carrier system in the absence of phonons is that of 
an Anderson insulator, i.e., zero. Phonons cause transi- 
tions between localized carrier states and, hence, allow 
conduction. Second, phonon-assisted carrier transitions 
meet the need to have carrier transitions between occu- 
pied states (which lie below the Fermi level) and unoccu- 
pied states (which lie above the Fermi level) whilst sat- 
isfying the demand that energy be conserved (Fig. p^ . 
This energy-conservation requirement also holds for the 
Hall effect. Now, as inelastic processes are being invoked 
one should ask the question: How can interference be- 
tween distinct hopping paths, necessary for the sensitiv- 
ity of the hopping rate to any quantal phase, arise? 

We will answer this question both in the context of 
longitudinal conductivity and the Hall effect. In the next 
subsubsection we will consider the longitudinal conduc- 
tivity. After that, in further sections, we will discuss 
what inelastic processes contribute to the Hall effect, and 
interference of those processes. 



3. Quantum interference in the presence of inelastic 
scattering. Longitudinal conductivity. 

In order to answer the question posed in the previ- 
ous section let us first formulate precisely what is meant 
by coherence and sensitivity to the quantal phase in the 
system at hand. The Hamiltonian of the system reads 



2. Hopping transport: compatibility of inelastic processes 
and quantum interference 

Despite the apparent simplicity of the foregoing anal- 
ysis of a triad of sites, the task of obtaining the linear 
dependence of the Hall resistivity on the magnetic field 
via Holstein's approach is a much more subtle matter. 
Furthermore, the issue of establishing quantal interfer- 
ence effects in this setting of inelastic hopping processes 
is equally subtle, so we shall now revisit this subject. 



FIG. 12. Schematic depiction of phonon-assisted inelastic 
hopping. Direct hopping (left) and indirect coherent hopping 
via an intermediate site (right) are shown. 



H — Hci + i?ci-ph + Hpii, (3.6) 

where i?ph is the phonon Hamiltonian. The exact 
quantum- mechanical eigenstates of i?ei, both localized 
and delocalized, are sensitive to applied magnetic field. 
(Here, we are, of course, concerned with localized states.) 
In the presence of a magnetic field, the exact eigenfunc- 
tions {r\'^j) can no longer be chosen to be real quanti- 
ties, and are thus characterized by an absolute value and 
a phase. It is often convenient to approximate the exact 
localized eigenstates in terms of the eigenstates {|(/>j)} of 
outer-shell electrons of isolated ions (where i enumerates 
the ions). We note that the {|0j)} are not, in general, 
orthogonal, although they are typically linearly indepen- 
dent and may, therefore, be taken as a basis. In terms 
of this basis a Hamiltonian describing carrier motion in 
the presence of the corresponding isolated ions located at 
positions {Rj} as well as a disordered potential has the 
form 
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El 



Vjk{4>k 



(3.7) 



where {ej} are random energies and Vjk are random 
transfer matrix elements. 

Having introduced the basis of locaUzed ionic states 
and the Hamiltonian for the disordered system 
of ions -ffcij we now examine in detail the effect on this 
system of an applied magnetic field. As shown by Hol- 
stein , in the presence of a magnetic field the ionic ba- 
sis states become the modified collection {\4'j)B}, 
being solutions of the Schrodinger equation 



1 



2ra* 



{Y\(j)j)B = e(r|0j)i3, 



(3.8) 



where U{y — Rj) is the potential of the ion located at 
position Rj, m* is the effective mass, and A is the mag- 
netic vector potential. The wave function {Y\4>j) b\b=o 
has the property that it is a function of r — Rj , and we 
would like to recover something like this property in the 
presence of the magnetic field. To do this we introduce 
a gauge transformation 



r K, 



(3.9a) 



T^j[T) = (e/?ic)A(Rj) • r = (e/?ic)(B x R^) • r, (3.9b) 

where we have chosen the gauge potential to be A(r) = 
H X r/2. It is straightforward to show that this trans- 
formation leads to an equation for (r|Kj) in which the 
vector potential (magnetic field) term contains coordi- 
nate in combination (r — Rj): 



^ (p + eH X (r-Rj)/2c)^ 



2ra 

C/(r-R,)] (r|0,)B 



(3.10) 



Hence, the wave function (r|Kj) is seen to have the sought 
property: 



{Y\nj)B = (r-Rj|Kj 



(3.11) 



It is convenient to employ the magnetic-field dependent 
ionic states {|(/'i)B} in the perturbative construction of 
the eigenstates of the system Hamiltonian (|3.7|) 



i/ei(B) = ^ \ct>,)B e,{B) {4>,\b + V^.fe(B) {4>k\B. 

(3.12) 

In what follows we shall omit the explicit dependence on 
B. 

We now construct approximations to the exact eigen- 
states of iJci(B) in terms of linear combinations of 
ionic states {\(l)i)B}- To do this we use renormalized 
(i.e. Brillouin-Wigner) perturbation theory in powers of 
Vkj/{Ej — Ek) (see, e.g., Ref. |67| , |68| ), thus obtaining 



E 



Ej — Ek 



VkhVhj 



{Ej - Ek){Ej - Eh) 



(3.13) 



where {Ej} are the exact energy eigenvalues. We now 
pause to remark that in the theory of hopping conductiv- 
ity in doped semiconductors ||5^,^,^,^ the parameter 
Vkj/{Ej — Ek) is indeed small, due to sizable distance 
between donors (or acceptors). In subsection devoted 
to application of hopping conductivity model to man- 
ganites we we will see that this parameter can also be 
rendered as small, because of the dependence of the ef- 
fective hopping amplitude on the core-spin misalignment. 
It is worth mentioning, however, that the sensitivity of 



I'^j) to phases arising from transformation Eq. (3.9a) is 



a general prope rty which does not rely on perturbation 



expansion 3.13 . It is also reasonable to assert that the 
essential dependence of on the magnetic field en- 

ters solely through such phase acquired by the local ba- 
sis wave fun ctions {{r\(j)j)} under the gauge transforma- 
tion in Eq. ( |3.9a ), because, at reasonable experimental 
strengths of the magnetic field much less than a quan- 
tum of flux treads a localized ionic wave function. 



Having constructed the states i.e., approximations to the exact localized states, we now use them to compute 
the square modulus of the matrix element Ujk of Eq. (3.4): 



|[/,fe|2 = |(v|/,|i7,i_ph|*, 



E E {(t^k\Hcl-ph\(f>j){(t)j\Hcl-ph\(t}k) 



VikVknVn 



i^k n^k,n^i 



{E, - Ek)HE^ - E„ 



- E E ('^i|-^ol-ph|0m)(0m|^^ol-ph|</'i 



VmjVjlVlm 



{E, ~ Ern)HE, ~ El) 



-fee. -I- 



(3.14) 



which features in the transition rate, Eq. (3.5). 
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Terms exhibited in Eq. (3.14) involve motion along 
paths that surround loops of nonzero area, i.e., the ma- 
trix elements of electron-phonon interaction and tra nsfe r 
amplitudes are taken between localized orbitals Eq. ( |3.§| ) 
of sites that form such loops. Note, however, that not 



all such terms are included in Eq. ( 3.14 ), but only those 
in which matrix elements of electron-phonon interaction 
enter the corresponding expressions in combination with 
their complex conjugated (i.e. time-reversed) counter- 
parts. Furthermore, the remaining matrix elements en- 
tering terms featured in Eq. ( |3.14 ), i.e., the overlap in- 
tegrals Vhk, involve motion along paths that surround 
loops of nonzero area . It is through such products of 
three overlap integrals VhkVkjVjh that the transition rate 
acquires its sensitivity to fluxes through loops of nonzero 
area. It is only such terms that lead to non-vanishing 
interference contribution to hopping probability which 
is sensitive to fluxes. Such sensitivity to phase results 
in the Aharonov-Bohm magnetoresistance effects in hop- 
ping conductivity (see, e.g., Ref. |f70[|). 

Let us discuss the physical meaning of interference 
terms featured in Eq. (3.14). Consider, for example. 



charge carrier hopp ing in a triad of sites. The relevant 
terms in Eq. ( 3.14 ), which involve loops, correspond to 
i = 1, j — 2, k — 2, n — 3 {in the first of featured terms), 
and m = 1, Z = 3 (in the second of featured terms). In 



terms of isolated orbitals the first term features in (3.14) 
corresponds to interference of two processes of charge car- 
rier tunneling from site 1 to site 2, a direct one and an 
indirect one (via site 3), with the following charge car- 
rier interaction with a phonon at site 2; the second term 
corresponds to carrier interaction with phonon at site 1 
with the following interference between direct and indi- 
rect tunneling paths from site 1 to site 2. (We note that 
the situation can, of course, be readily generalized to 
the case of changes in more than one phonon occupation 
number.) 

We now notice that by contrast with terms that are 
featured in ( |3.14 ), terms that have been omitted there, 
like 



siUjki"^ = y^(0i|gci-phi'/'j)('/'j|gci-phi0fc 



Ei — Ek' 



(3.15) 



involve motion along paths that surround loops of 
nonzero area, but do not lead to non-vanishing interfer- 
ence contribution to hopping probability. The reason for 
absence of such contributions is that they do not involve 
combinations of matrix elements of electron-phonon in- 
teraction with their time- reversed counterparts. To 
see why such combinations are important, consider the 
electron-phonon interaction Hamiltonian 



l-ph 



(3.16) 



where q is the phonon wavevector. Evaluation of non- 
vanishing interference contribution, therefore, includes 
sums over all phonon wavevectors. Terms that exhibit 
two matrix elements of electron-phonon interaction that 
are not complex conjugated to each other turn out to 
be oscillating functions of q, and vanish upon evaluation 
of the sum over q. This simple observation allows us 
to find all important interference terms in the hopping 
probability (3.14 without resorting to explicit evalua- 
tion of {(f)i\Hc\--p\-i\cj)'j) as it was done, e.g., in j2^ in 
consideration of the Hall effect and in |Q in considera- 
tion of the Aharonov-Bohn hopping magnetoresistance. 
We note that if the sum of two terms featured in (3.14) 
is equal to zero, this consideration easily allows us to 
find the appropriate next terms of expansion of hopping 
probability in powers of Vij /{Ei — Ej). 

We are now in a position to discuss what is meant by 
phase coherence and sensitivity to quantal phases in the 
present setting of hopping conduction. In the absence of 
inelastic agents, products of, e.g., three overlap integrals 
VhkVkjVjh (where sites h, k and j form a nondegener- 
ate triangle, with non-zero flux threading this triangle 
in the presence of magnetic field), and carrier energies 



Ej (see the remark 



|) are sensitive to the Aharonov- 
As we can see from Eq 



Bohm quantal phase. As we can see from Eq. (3.14), 
interference of two amplitudes of quantal transition be- 
tween states h and k occurs (and these amplitudes are 
coherent) even if these transitions are due to inelastic 
agents. For such interference to occur, inelastic agents 
which determine one of the amplitudes of the transition 
must be the same as inelastic agents which determine an- 
other amplitude of the transition (i.e. phonon frequencies 
and wavevectors are equal), so that the square modulus of 
matrix element of, e.g., the electron-phonon interaction, 
enters the probability. As we have mentioned above, the 
presence of square modulus of the electron-phonon inter- 
action matrix element in the probability corresponds to 
phonon factors in the two interfering amplitudes (e.g., 
in the amplitude of a direct process and in the complex 
conjugate of the amplitude of an indirect process) being 
related to each other by the time-reversal. 

Let us now discuss what is meant by the phase breaking 
or decoherence in the context of the hopping conductiv- 
ity regime. The eigenstates and their energies Ej 
are determined by 7Jei-ph, i-e-, in the absence of phonons. 
When phonons interact with charge carriers, one can con- 
sider them as a reservoir which leads to randomization 
of carrier states. The time scale of such randomization 
is given by the inverse rate of carrier transitions between 
different eig enst ates of Hd-ph caused by electron-phonon 
interaction ( |3.5| ). It is this randomization that is meant 
by phase breaki ng ( or decoherence). However, the carrier 
transition rate ( |3.5[ ) between different states itself, which 
measures the rate of decoherence (being determined by 
eigenfunctions of these states) carries information about 
quant al phases that arise in the context of Eqs. (3.13) and 
( |3.14 ). Furthermore, the electric current arises during 
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the act of hopping, i.e. when randomization has not yet 
occurred. Therefore, the answer to the question, posed in 
the second paragraph of the present subsection, concern- 
ing interference in the presence of inelasticaUy-assisted 
hopping is as follows: The (steady-state) hopping current 
is generated during the process of inelastic scattering, 
whilst decoherence arises only after this scattering has 
occurred; thus, decoherence effects do not preclude sen- 
sitivity of hopping conduction to quantal phases. More- 
over, the question of interference of consequent hopping 
events is meaningless in the context of hopping conduc- 
tivity, because only amplitudes of hopping between the 
same initial and the same final states can interfere. Thus, 
there is a significant difference between interference ef- 
fects in hopping conductivity and in diffusive mesoscopic 
transport. In diffusive transport, the whole diffusive tra- 
jectory, with all consequent scattering events, determines 
the current via the diffusivity, and electronic coherence 
during consequent scattering events is important for in- 
terference effects. In hopping transport, incoherence of 
consequent hopping events does not contradict interfer- 
ence of quantum-mechanical amplitudes that determine 
a single charge carrier hop and the hopping current itself. 



4-. Quantum interference in the presence of inelastic 
scattering; processes leading to the Hall effect 



For localized carriers, the interference processes that 
lead to the Hall conductivity are even more peculiar than 
those leading to the sensitivity of the longitudinal con- 
ductivity to the Aharonov-Bohm phase. In particular, 
the occurrence of the Hall effect requires phonon-assisted 
hoppings between exact initial and final carrier states, 
via an exact intermediate state. (It is not sufficient to 
include as intermediate states the virtual ionic orbitals 
that enter via Eq. (|1|.) The reason for this difference 
between the hopping Hall conductance and the sensitiv- 
ity to magnetic flux of the hopping magnetoconductance 
arises because of the necessity to extract a dependence of 
the Hall conductance that is linear in (or, more generally, 
an odd power of) the magnetic field . 

Let us now explore the collection of processes that con- 
tribute to the Hall conductance. These processes, first 
identified by Holstein, can be determined by making use 
of the odd (i.e. dissipative) character of transition rates 
under the time-reversal operation: t — > —t. As discussed 
in Sec. HI C and HI D, the hopping current is determined, 
via the conductivities of the resistive network, by the 
rate of transitions between sites. The Hall current, as 
with any current, is odd under time-reversal. So, too, is 



the transition rate. (To appreciate this, consider an el- 
ementary relaxation process for some quantity Q, which 
evolves according to the rate equation 



dQ 
dt 



Q 



(3.17) 



From the consistency of equation one immediately sees 
that it is formally appropriate to designate the relax- 
ation rate r as being odd under the time reversal.) Now 
consider the transition probability per unit time between 
the exact single-carrier states i and /, viz. Wfi which, 
according to Fermi's Golden Rule, is given by 



2ti 



A},\H{E,-Ef), 



(3.18) 



where Afi is the sum of the transition amplitudes for all 
coherent (i.e. interfering) processes connecting i and /, 
and the 5 function imposes energy conservation between 
the initial and final states, the energies Ei and Ej of these 
states being full energies. (That is, they include not only 
carrier but also phonon energies) . The (5-function, which 
has the complex representation 



1 . i 

S{E) = Im— lim , 

TT s^+o E — is 



(3.19) 



is an odd quantity with regard to time-reversal, in the 
sense that under the transformation t —t, sign of imag- 
inary part s in the denominator changes and, thus, so 
does the sign of the S function. By contrast, the quantity 
is even (i.e. non-dissipative) under time-reversal. 
We note, in passing, that the precision of the energy con- 
servation is not what is essential (from the point of view 
of time-reversal properties). For example, the imaginary 
part of a Lorentzian function. 



Im — 



1 



TT E-iT TT E^ +T^'' 



(3.20) 



reveals that this function, too, inherits the oddncss of the 
rate F. 

En route to exploring the processes that contribute to 
the Hall conductance, let us now consider a simple case 
in which Afi has just two contributions: 



A ... — /(dir I 4ind 



(3.21) 



where A'jf is the amplitude for the direct path and A™f 
is the amplitude for an indirect path (i.e. a path via an in- 
termediate state). These two amplitudes can be written 
in the form 



^dir _ ^0,dir 
yiind 



A 



0,ind 



fi 



(3.22a) 
(3.22b) 
(3.22c) 



where Aj:^ " and Aj:\ " are the zero-magnetic-field ampli- 
tudes, and and $2 are phases arising in the presence of 
magnetic field for direct and indirect paths, correspond- 
ingly, and 4> is the difference of phases between direct 
and indirect paths induced by magnetic field, i.e., the 
Aharonov-Bohm phase. (We neglect changes magnetic- 



field-induced changes in the magnitudes of A^jrf" and 
^o.dii -J p^^. transition probability we then have 
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HI 2 



0,ind|2 



1^' 

him A 



Re A^. 



A^j ™'^ cos( 



O.dir * ^O.ind 



(3.23) 



We now observe that the term containing sin is the only 
contribution to the probabihty that is odd with re- 

spect to the transformation cj) —cj) (i.e. with respect to 
magnetic- field reversal), this oddness being a necessary 
property of the Hall conductance. Thus, in a computa- 
tion of the Hall conductance, only the imaginary part of 
the quantity corresponding to A^jf^^ contributes 
and, therefore, one has to consider those indirect pro- 
cesses for which the zero-magnetic-field amplitude has 
a component out-of-phase with the zero-magnetic-field 
amplitude of the direct process. (Said another way, one 
must consider contributions to A^rf" * A^r'^'^'^ that are odd 
with respect to time reversal.) Such contributions do not 
appear if, as in the case of the longitudinal hopping con- 
ductivity in Sec. HID 3, one considers one-phonon pro- 
cesses They do, however, appear if one considers, 
e.g., two-phonon processes described by combinations of 
amplitudes obeying the following property: A'jf" con- 
tains an even (odd) number of complex energy denomi- 
nators when A^^'j'" contains an odd (even) number. Then 
these denominators give rise to an additional energy- 
conserving ^- function [ fz^ , a q uantit y that is odd re- 
spect to time reversal [cf. Eq. ( 3.19| )], and yields con- 
tributions to the probability that behave suitably under 
time-reversal. 
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FIG. 13. Pairs of interfering amplitudes that result in the 
Hall effect. Rows (a) and (b) each show the interference of 
two-stage processes. Rows (c) to (f) each show the interfer- 
ence of a one-stage and a three-stage process. Right column: 
indirect hopping processes; left column: direct hopping pro- 
cesses. Lines with arrows correspond to carrier propagators; 
their intersections with wavy lines correspond to matrix ele- 
ments of carrier-phonon interactions. 

Having ascertained the structure of the amplitudes 
that give rise to the Hall effect, we now select and ex- 
amine the dominant contributing processes (i.e. those in- 
volving the smallest possible number of electron-phonon 
interactions). These involve two-phonon transitions and, 
as shown by Holstein their elements (i.e. direct and 
indirect processes) can be visualized as follows: 
(i) Both transitions are two-stage processes (i.e. involve 
two transitions between exact carrier states. For exam- 
ple, the indirect and indirect transitions respectively be- 
ing 



(fc,7V{,iV^), 
{k,Ni,N^), 



(3.24a) 
(3.24b) 



where A^i, N[, N2 and N2 are the occupation numbers 
of phonon modes 1 and 2, N, and N' differ by unity, i, 
j and k label initial, intermediate and final sites (as well 
as carrier states localized on these sites), 
(ii) The direct transition is a one-stage process, e.g.. 



{i,m)^{k,N[), 



(3.25a) 



and the indirect transition is a three-stage process, e.g.. 



ii,Ni,N2)~^ ii,N[,N2) 

^ {j,N[,N^)^ik,N[,N2). 



(3.25b) 



(The whole set of processes includes those which re- 
sult from alterations of the sequences of various subpro- 
cesses). One can observe, that phonon modes whose pop- 
ulation is changed in processes (i) or (ii), interact with 
charge carrier on both direct and indirectpath, and, as 
we shall see, these paths can interfere [Q. Complete 
set of processes that lead to the Hall effect is shown on 
Fig. |ll. 



In terms of phonon-assisted transitions between exact carrier eigenstates the two-stage processes are described 
by the amplitudes 

.dir,2 _ 

^k,Ni,N2;j,N[.N.^ ~ 
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(3.26a) 



Ej - Ek 



E , - E, 



<ind,2 



E,, - Ej + {N[ ~ Ni)hLOi + ih-f E, - Ej + {N^ - N2)huj2 + ih"f 

where ^^'^^ N2-j n' n' ^^'^ ^^ti^ N2-j n' n' respectively, the ampHtudes of the two-stage-dkect and two-stage- 
indirect processes. The interfering amphtudes for the one-stage process A*^"'^ and (an example of) a three-stage 
process A"^'^'^ are given in terms of the exact carrier eigenstates 

A"""^^ = {^k,Nl\H,l^ph\■^^,N[), (3.27a) 

,ind,3 ^ (^fc, iV2|gcl-ph|^,, iV^|gcl-ph|^j-, A^2) iVl|gcl-ph|^., Nj) 

{E, - Ejm - Ej + {N[ - Ni)nw^ + ih-i) ■ ^ ■ > 



Let us briefly discuss the energy conservation (i.e. 5- 
function) structure and time-reversal properties of the 
probabilities associated with the amplitudes given in 
Eqs. (|j26|), (|^), ( F27al) an d ( |3l7b| ). The explicit 
(5-function in the formula ( 3.18 ) for Wft constrains the 
energies of the initial and final states. One further 6- 
function arises when we insert explicit ex pressions for 
the interfe ring pairs o f amplitud es, E qs. ( 3. 26a . 3. 26b ) 



and Eqs. ( |3.274|3.27h| ), into Eq. ( p^ . This second 5- 



function characterizes energy conservation between ini- 
tial and intermediate (or intermediate and final states). 
As we mentioned earlier, time reversal symmetry can also 
be satisfied if we use, e.g., imaginary parts of Lorentzian 
functions instead of (S-functions, thus taking into account 
approximate energy conservation in transitions between 
broadened states of the system. Here, for brevity, we use 
the term (5-function when discussing time-reversal prop- 
erties of transition rates. 

Let us follow how these 5-functions (and, hence, the 
requisite odd character under time reversal) emerge. For 
two-stage processes, a direct transition has to contain one 
intermediate state, which has to be virtual state, and. 



therefore Eq. ( 3.26a ) contains real energy denominators. 
Thus, in two stage processes, a i5-function additional to 
one giving energy conservation between initial and fi- 
nal state, is to arise from indirect transition amplitudes. 



Such a (5-function indeed arises (in each of the contribut- 
ing terms) from complex denominators of Eq. (3.26b), 
upon summation over all possible phonon modes. In in- 
terference involving three-stage processes, the a nalytical 
expression for three-stage amplitude Eq. ( 3.27b| ) is char- 
acterized by two energy denominators. Upon summation 
over all possible phonon modes, one of these denomina- 
tors lea ds to a iJ-function. Another energy denominator 
in Eq. (3.27b) corresponds to a virtual transition. As 
follows from Eq. ( 3.23 ), in the presence of the Aharonov- 
Bohm phase 4> picked up by carriers moving around three 
sites, interference contributions to Wj f det ermined by 
Eqs. ( |3.264|3.26b| ) and by Eqs. ( ^.27at|3.27b| ) all contain 
two (5-functions, and, as required by time-reversal sym- 
metry properties, are proportional to sn\<f>. 

In order to see the physical meaning of interference 
between amplitudes of direct and indirect two-stage pro- 
cesses or interference between amplitudes of one- and 
three-stage proces se s, it is in struct ive to write down the 
amphtudes ( |3.26a| ), ( ^.26bD , ( |3.27a| ) and ( |3.27bD in terms 
of the ionic orbitals \(j)i) and the transfer integrals Vy , 
by using Eq. (3.13). (The energy denominators are cor- 



rected compared to those that can be found in the orig- 
inal Holstein paper |2^). This will also allow us to for- 
mulate conditions necessary for coherence of two-phonon 
processes, relevant for the Hall effect. 



In terms of and Vij, the two-stage processes are described by the amplitudes 

.dir,2 

^k,Ni,N2\j,N[,N'^ ~ 

((/.fc,A^l|gel-ph|'/'fc,A^{)('/'»,A^2|gcl-ph|4,^2)^fc»-('/'fc,^2|gcl~ph|0fc,^ 

{E,~Ek){E,^Ek) 

, (0fc,7V2|i?el-ph|0fe,Af^)(0z, Afl|i?ol~ph|(/'«,iV{)l^zfe - ((/.fe,iVi|i/cl-ph|0fc,iV{)((/.„iV2|i/ol^ph|0z,Af^)l^.fc 



{Ei — Ej){Ek — Ei) 



(3.28a) 



A 



id. 2 



k,Ni,N2\j,N[,N'^ 



17 



{E, - E, + {N{ - Ni)hLOi + ih-i){Ek - Ej) 



E, - Ej + {N!^ ~ N2)huj2 + ih-f 



(3.28b) 



The interfering amplitudes for the one-stage pr ocess A '^"'^ and a n example of a three-stage process A'"'^''^ which, in 
terms of exact carrier states, are given by Eqs. (3.27a and (3.27b) in terms of ionic orbitals read: 



^dir,l ^ 
^ind,3 



ik,Nl\H,l^pl,\^k,Ni)Vk^ + (0„iVi|i/el-ph|0»,A^()^fcz 



+ 



Ei — Ek 

{Ej - Ek){E, - E,){E, - Ej){E, - E, + {N[ - Ni)huji + th^) 



(3.29a) 
(3.29b) 



First term in Eq. ( p.29a| ) does not contribute to inter- 
ference of one- and three-stage processes, because the 
matrix element {(j)k, Ni\Hei-pii\(j)k, N{) does not corre - 



spond to any time-reversed counterpart in Eq. (3.29b). 
We are now in a position to discuss the physical meaning 
of processes that contribute to the Hall effect in terms 
of local orbitals. In the interference of two stage pro- 
cesses, both amplitudes, direct and indirect, correspond 
to interaction with a phonon mode {Ni,uji) at site k 
(initial state |'/>fc), tunneling to site i (final state \(f>i) ) 
directly or via intermediate site j, and interacting with 
a phonon mode {N2,uj2) (at site i) that is distinct from 
one participating in a process occurring at site k. In 
the interference of one- and three-stage processes, the di- 
rect one-stage process (which is, strictly speaking, can be 
called one-stage only in terms of exact carrier states) can 
include two possibilities: (i) Interaction with a phonon 
mode (NijUJi) at initial site and tunneling to the final 
site; (ii) tunneling to the final site and interaction with 
a phonon mode {Ni,uJi) at the final site. Then the three 
stage process have to include the following stages: In 
the case (i) there is interaction with the phonon mode 
{Ni,uJi) at initial site, tunneling to intermediate site j, 
emission (absorption) and reabsorption (reemission) of a 
phonon mode (^^2, ^2) at site j, and tunneling to the final 
site. In the case (ii) charge carrier tunnels to intermedi- 
ate site j, where emission (absorption) and reabsorption 
(reemission) of a phonon mode (A^2, ^^^2) occurs, and then 
tunnels to the final site. We therefore see that, in terms 
of local orbitals, processes contributing to the Hall effect 
are characterized by interference of amplitudes in which 
phonon modes changing their occupation numbers are 
represented by time reversed counterparts in A'^" and 
A™'^, respectively. 

Having described the inelastic processes leading to the 
Hall effect, we are now in the position to generalize 
conditions for occurrence of interference, and to formu- 
late these condi tions for two-phonon pro cesses. It fol- 
lows from Eqs. (|3.26b| , |3.26a| , |3.27a| , |3.27b|) that electron- 
phonon interaction results in coherence of transfer ampli- 
tudes in two cases: (i) If direct transition and transition 
via an intermediate state both occur as two-phonon pro- 
cesses , with two phonon modes changing their occu- 
pation numbers in the course of both these transitions, 
then interference exists when the two phonons leading to 



the direct transition are the same as two phonons leading 
to a transition via intermediate site; (ii) If the direct tran- 
sition occurs as one-stage process assisted by a phonon 
mode, this phonon mode changes its occupation number, 
and the transition involving an intermediate site occurs 
as a three stage process. Then, the condition is that 
one of the phonon modes assisting three-stage process is 
the same as that in the direct transition, while another 
phonon mode which assists the indirect transition does 
not change its occupation number. 

Broadly speaking, both these cases lead to conditions 
that inelastic modes that change their state in the course 
of electronic hop are the same in the two interfering hop- 
ping amplitudes. Only in this case an interference exists 
even between amplitudes of inelastic processes. We note 
that this point has been also recently revisited by Entin- 
Wohlman et al. l75|l. 



5. Ordinary Hall effect: Local conductivities. Remarks on 
averaging over triads 

These interference contributions will result in the or- 
dinary Hall effect, with the Hall conductivity in a triad 
of ions CToH given by 



o-QH = G{ej} sin (B • Q/(f>o) 



(3.30) 



where 4>o is the (electromagnetic) flux quantum, Q is 
the (oriented, real space) area of the triad, and 
are the energies of the three single-particle eigenstates, 
which are invariant under reversal of the AB flux. The 
explicit expression for G can be found by substituting 
Eqs. (gj28|), ( |3^28b| ), ( p9al ) and ( |l29b| ) into Eqs. ( |3l^ ) 
and ( |3.3| ) . (Note that in generic case G also depends on 
the populations of these states, which themselves may 
depend on particle-particle correlations.) 

In j2^, Holstein mainly addressed the issue of the 
Hall effect in hopping conductors in the presence of an 
a.c. electric field. Compared to the d.c. Hall effect, the 
a.c. problem is simplified because the principal contribu- 
tion to the Hall conductance comes from those spatially 
isolated configurations of sites for which the population 
relaxation time is on the order of the inverse frequency 
uj of the current (i.e. ujti- ~ 1). In this a.c. case, there is no 
need to address the (highly nontrivial) issue of how these 
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sources of the HaU effect (i.e. configurations of sites) are 
combined into a conducting network that is connected 
to the Hall contacts. For the d.c. Hall effect, on the 
other hand, this issue of the structure of the conduct- 
ing network must be faced, and it becomes necessary to 
understand which triads are the most effective in con- 
tributing to the Hall effect, how to average over triads, 
and what quantity should be averaged. [If the quantity 
to be averaged should be the conductivity (resistivity) 
then one should first compute the local conductivities 
(resistivities) of triads and then obtain the macroscopic 
conductivities (resistivities) by averaging.] This issue of 
averaging over all triads and conducting network struc- 
ture in disordered systems remains controversial ||7j,^. 

However, in the present Paper we are concerned not 
with the ordinary Hall effect in a system with localized 
carriers but the anomalous one. For the Anomalous Hall 
effect, the Aharonov-Bohm phase does not play a fun- 
damental role, a very weak magnetic field being applied 
solely for the purpose of inducing a macroscopic magneti- 
zation of a ferromagnetic medium. Rather, for the AHE 
in system comprising magnetically disordered core-spins 
on Mn sites visited by hopping charge carriers, it is a cer- 
tain type of spin quantal phase [|6 47 4^ that manifests 
itself. We now turn the origin and meaning of this spin 
quantal phase. 



E. The quantal Pancharatnam phase 




1 

O 2 



FIG. 14. (a) Triad of Mn ions having distinct core spin 
orientations, (b) Sphere of possible core spin orientations, 
showing the specific orientations of the spins for this triad. 
The orientations form the vertices of spherical triangle. The 
area of this triangle determines the quantal Pancharatnam 
phase. 

To understand the nature of the spin quantal phases, 
let us begin by examining the single-particle quantum 
mechanics of a carrier hole added to a triad of Mn^^ ions. 
We regard the spin-3/2 core spins of the Mn ions as large 
enough to be treated classically, so that one can assign a 
definite direction to each of them. Thus, a generic config- 
uration of core spins is characterized by the triad of unit 
vectors {ni,n2,n3}, respectively located at the triad of 
sites {Ri,R2,R3}, as depicted in Fig. [l^. Let us now 
consider the transfer of the hole carrier between ions. In 



a, 3 



the double exchange model, such transfer is described by 
the Hamiltonian 

+J \a,(t)j) Uj ■ cTa^p {PAjl (3.31) 



where cr is the Pauli spin operator describing the spin 
of the hole carrier, J is the Hund Rules coupling energy 
which, for Mn"^^ ions, is on the order of several eV, and is 
much larger than the orbital transfer integrals Vjfe, and 
|q;,(/)j) is the outer-shell atomic carrier state at site j, 
the index a labeling the spin-projection on to the z axis 
(i.e. |q!)|q=+ areeigenstatesof (Tz). In general, the Hamil- 
tonian ( 3.31 ) results in spin polarization of electrons at 
arbitrary ratio of J and Vju- However, two limiting cases 
are of interest. At J <C \Vjk\, metallic ferromagnetism is 
usually treated perturbatively in terms of an electron gas 
in metals, resulting in the RKKY interaction. In this case 
charge carriers with any spin projection are taken into 
account. The opposite case J ^ \Vjk \ is a purely double- 
exchange model. This is the case that is relevant to man- 
ganites. By Hund's Rules (appropriate for J ^ l^fel), 
there is, at each site, effectively a single quantum state 



available to the hole. This state In, 



is the one in 



which the carrier spin projection opposes the core spin 
direction. The orbital \(t)j) that characterizes this state 
is one of the orbitals of an isolated Mn ion, centered on 
the site j; we choose to omit the remaining orbitals so 
as to simplify the discussion. As for the other spin state 
(as well as any other orbitals), we regard them as being 
inaccessible on energetic grounds. Postponing until later 
any effects of spin-orbit interactions, we assume that the 
transfer of carriers (being, as it is, effected by either the 
"kinetic energy" by the second term in Eq. ( |3.31 ), or by 
the electron-phonon interaction) has no effect on the spin 
of the carriers. The hopping amplitudes between ionic 
states are thus given by (n^, 0fe|r|nj-, c/f^), where T corre 



sponds either to the transfer operator V of Eq. ( 3.7 ) or to 
the electron phonon interaction iJei-ph- These hopping 
amplitudes depend explicitly on the relative orientation 
of the core spins, a nd rij. In particular, by projecting 
the Hamiltonian ( 3.31 ) on to the physically relevant low- 



we arrive 



energy subspace spanned by the states |n 
at the the double-exchange Hamiltonian projected on to 
the low-energy subspace, -ff^E ; which takes the form 



(3.32a) 



V,k [ cos y cos ^ + e'-"-' sin y sin ^ ) , (3.32b) 
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where V^^ are the effective transfer amphtudes, jkj = 
7fc — 7j" and, respectively, Oj and 9k are the azimuthal 
angles and jj and jk are the polar angles of the semi- 
classical spin directions rij and Uk- Provided we choose, 
e.g., rij II Bz and || e^, (where {ej:,ey,ez} are the 
Cartesian basis vectors) this effective transfer amplitude 
reduces to the Anderson-Hasegawa form: 



— = tan 
2 



ni • (n2 X ng) 



KtH = F,feCos0/2, 



(3.33) 



where 6 is the angle between rij and n^. However, and 
this is central to our discussion, the effective transfer am- 
plitude is, in general, a complex quant ity ch aract erize d 
by its amplitude and phase. From Eqs. ( 3.32b and ( 3.33 ), 
it is indeed apparent that if core spins are co-aligned, the 
effective hopping amplitude is maximal, while if the core 
spins on two ions are opposite, hopping between such ions 
does not occur. 

We now follow the line of argument applied in 
Sec. Ill D 3| , and co nstruct the exact eigenstates of 
Hamilitonian ( 3.32a ), via Eq. ( 3.13 ) modified to ac- 
count for spin. Hence, we can build matrix elements 
of the electron-phonon interaction between exact local- 
ized states, taking into account the effect of the core spin 
orientations {rij}. We now note that this expansion for 
the state gives rise, in the hopping probability |J7jip (and 
hence in the hopping rate Wji), to terms containing prod- 
ucts of matrix elements such as transfer integrals V^^ 
and electron-phonon interactions Ukj- Amongst these 
terms are ones containing matrix elements associated 
with paths around closed loops and incorporating the ef- 
fects of interference between distinct carrier paths. The 
simplest example involves the product V^^V^^Vj'^^ , as- 
sociated with the path i ^ j ^ k ^ i. 

In the presence of the constraints set by the core spin 
orientations, the transfer of carriers discussed in the pre- 
vious paragraph is subject to a striking quantal effect. 
To see this, consider the products of matrix elements 
^ik^kf^j^- Explicitly, such products have the form 

T/eff T/eff T/®^ 
^ ik kj *^ ji 

= (ni,0j|iJDE|nfc,0fc) 

X (tifc , 0fc liJoE I rij , Hde \ , (f>i) 

= (n,| (g) {(t)i\HijE\(j)k) ® Irifc) 



x(nfe| (g) (0fc|i?DEl<^j) ^ 



X n,- 



n, 



(X (n,|nfc)(nfe|n,)(n,|n,) - TrP^ P„ (3.34) 

where the operators Pj = {l + cr- nj)/2 are projectors (in 
spin space) on to the spin states aligned with the local 
core spin orientations n^. From this last expression, in 
terms of projectors, it is straightforward to establish that 

TrPfcP.P, (3.35) 
= (l ni • n2 + n2 • na -t- na • rii) + i(ni • (n2 x ng)). 

Hence, we arrive at the quantal phase f2, the phase of the 
complex quantity TrP^ Pj Pi, which is given by 



H- ni • n2 + n2 • na -I- ns • ni 



(3.36) 



In the context of the physical quantity from which 
the computation of the quantal phase f2 emerged, viz., 
the perturbative evaluation of the hopping rate Wkj, 
Eq. (3.18), this phase modulates the interference between 
hopping processes that progress from one site j to an- 
other site k, either directly or indirectly, via a third site. 



Formula Eq. (3.36) indicates that the phase 11 has a 



geometric interpretation as the (oriented) solid angle of 
the geodesic triangle having vertices at {ni,nj,nk} on 
the unit sphere. It is the quantal analog of the classical 
optical phase discovered in the context of polarized light 
by Pancharatnam j2^.2^. In that setting, what Pan- 
charatnam showed is that a under a sequence of changes 
of the polarization state of light that return the light to 
its original polarization state there arises a phase shift 
(i.e. a phase anholonomy) determined by the geometry 
of the sequence of changes. If the sequence of polariza- 
tion states is represented by a sequence of points on the 
Poincare sphere (a certain parameterization of light po- 
larizations) then n is given by the area of the geodesic 
polygon on this sphere the vertices of which are corre- 
spond to these polarization states. 

In the double-exchange electronic analog of Pancharat- 
nam's phase, the transporting of an outer-shell carrier to 
an ion with a differently oriented core spin via a spin- 
independent process is characterized by a matrix element 
that can be interpreted as a connection. For processes 
visiting a closed sequence of sites, this connection yields 
a quantal phase Q, viz., the phase shift of the return- 
ing spin state in terms of the sequence of sites visited. 
This quantal version of Pancharatnam's phase is given 
by half the area of the geodesic polygon (on the unit 
sphere of core spin orientations) the vertices of which are 
the core spin orientations of the sites visited. Although 
the phase f2 emerged from considerations of interference 
between processes involving a triad of sites, such phases 
are more general and would, in fact, emerge for arbi- 
trary processes. We remark that, in contrast to Berry's 
adiabatic phase | |50[ |, the phenomenon described here is 
associated with sudden changes in the carrier-spin state, 
and need not be slow. 



F. The Anomalous Hall effect in hopping regime 



We now turn from the OHE in a spinless triad to the 
AH E in a triad of magnetic sites. Like the OHE given by 
Eq. 3.30| , this AHE results from two-phonon processes, 
but is due to the Pancharatnam phase instead of the AB 
phase. (At this stage, we have not yet included the ef- 
fects of the spin-orbit interaction.) Mutatis mutandis, we 
arrive at the AH conductivity ctah, given by 
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1 ^13 ^32 ^21 . ^ I 

0-AH = C^{£il COS — COS — COS — sm — , (3.37) 

where cosOjk = rij ■ n^, the factors cos{9jk/2) are 
Anderson- Hasegawa factors, and {ej} are the energies 
of the three on-site single-particle eigenstates that are 
consistent with Hund's Rules, these energies depending 
on {rij • nfc}i<j<fe<3 and cos(0/2). Note that G is even 
under the reversal of the Pancharatnam flux —^ — il, 
and a AH is odd under it. 




FIG. 15. Two configurations [(a) and (b)] of core spins on 
a triad of Mn ions, differing by the interchange of the spins 
on sites 2 and 3. These distinct spin configurations give rise 
to opposite Pancharatnam phases. 

We have shown that, for a triad with given set of core- 
spin orientations, an AHE arises from the quantal Pan- 
charatnam flux. However, there is a significant difference 
between this AHE and the OHE. In the former (non- 
magnetic) case, a uniform applied magnetic field leads 
to a net macroscopic OHE, even though contributions of 
triads may cancel one another |7^. In the latter case 
(with its magnetic sites, Pancharatnam flux, but spin- 
orbit interactions not yet included), even the presence of 
a macroscopic magnetization of the core spins is insuf- 
ficient to cause a macroscopic Hall current. The reason 
for this is that in obtaining the macroscopic AH current 
from Eq. ( 3.37 ) we must average over the configurations 
of the core spins. In the absence of spin-orbit interac- 
tions, the distribution of these configurations, although 
favoring a preferred direction (i.e., the magnetization di- 
rection m = M/M), is invariant under a reflection of all 
core-spin vectors in any plane containing the magneti- 
zation, and, therefore, there is no preferred Pancharat- 
nam flux. For example, two spin configurations shown in 
Fig. ^ have the same magnetization but opposite signs 
of the Pancharatnam flux. This fact, coupled with the 
fact that {sj} are also invariant under such reflections, 
guarantees that the macroscopic AH current will average 
to zero. (We do, however, expect significant AH current 
noise in the ferromagnet-paramagnet transition regime, 
owing to the fluctuations of the Pontryagin charge ||7^] 
of the triads of core spins (which we shall define shortly) 
and, hence, elementary Pancharatnam fluxes.) 

In order to capture the AHE in double-exchange ma- 
terials such as manganites, we must consider not only 
the Pancharatnam phase but also some agent capable of 
lifting the reflection invariance of the energies {sj} and 
the distribution of core-spin configurations and, hence, 
capable of inducing sensitivity to the sign of the Pan- 



charatnam flux. Such an agent is provided by spin-orbit 
interactions. We now discuss the effect of these interac- 
tions on the motion of charge carrier in a triad. 



G. Spin-orbit interactions in a triad 

The most general form of spin-orbit interaction is given 
by the spin-orbit Hamiltonian 



Hso = ap • (cr X VU) , 



(3.38) 



where the potential U includes ionic and impurity po- 
tentials, a is the spin-orbit interaction constant, p is 
the electron momentum, and cr are the Pauli operators. 
This spin-orbit interaction results in an effective SU(2) 
gauge potential Ago = am{(T x V?7) |Q, where m is the 
relevant mass of the carrier. This gauge potential pro- 
vides an additional source of quantal phase. For a given 
core-spin configuration, the spin-orbit interaction favors 
one sense of carrier-circulation around the triad over the 
other, and thus favors one sign of Pancharatnam phase 
over the other. 

Let us consider the consequences for the energy spec- 
trum of the triad {sj} that arise due to spin-orbit interac- 
tions. This interaction generates a dependence of {sj} on 
the three vector-products Nj-fe = rij x rife which, together 
with the magnetization direction m, yield a preferred 
value for the triad Pontryagin charge qp [= rii • (n2 x n^)] 
and, hence, a preferred Pancharatnam flux. 

It is straightforward to find corrections, due to the SOI, 
of hole eigenenergies if the on-site energies of the holes 
are nondegenerate. Then the sensitivity of {ej} to vector- 
products Njfe = rij X rife first enters at third order (in the 
transfer matrix elements): 



-eh){s,-ek), (3.39) 



where Tjk = Pj Vjk Pk are the transfer amplitudes, Vjk 
are the hopping matrix elements, and Tr denotes a trace 
in spin space. (For degenerate on-site hole energies one 
should obtain the splitting of these energies due to trans- 
fer in the absence of SOI, and then include SOI at the 
final step, arriving at the result to be given below.) The 
hopping matrix elements are sensitive to the SOI quantal 
phase, and can be written in the form Vjk = Vjk^ Ljk, 
where Ljk = (l -I- icr • gjk), Vjk^ is an orbital factor, and 
gj-fe (oc ttso) is an appropriate vector that describes the 
average SOI for the transition j — > fc in a triad of sites i, 
j and k. Then, e.g., the first-order (in a) shifts in the e's 
are given by 

Ssj (X Tr Ti3 T32 T21 = 4 Re Tr Pi L13 P3 L32 P2 L21 

= 2 (Ni3 • gi3 + N32 • g32 + N21 • g2i) - N • g, (3.40) 

where N = N13 + N32 + N21, and g = gi3 + g32 + g2i- 
If the potential U in the SOI is a superposition of 
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spherically-symmetric ionic potentials in a the triad of 
sites then the vectors gjk have a transparent geometrical 
meaning: 



Q 



1 

2 



(R. 



Rh) X {Rk - Rh) , 



(3.41) 
(3.42) 



i.e., they are proportional to the area |Q| of the triangle 
whose vertices are the sites Rj, R^ and R/j. Then the 
SOI-generated shift in the carrier eigenenergies has the 
Dzyaloshinski-Moriya form . 



H. Elementary Hall conductivity in a triad 

There are two contributions to the AHE which result 
from the SOI-generated shift in the carrier eigenener- 
gies. The first contribution is due to the dependence 
of the probability of hopping around the triad on {sj} 
for a given spin configuration. By incorporating the 
shifts ( |3.40 ), together with the Pancharatnam phase, we 
arrive at the elementary AH conductivity 



crj^y = ni • (n2 X na) Se^ dG/dsj. 



(3.43) 



As discussed above, Eq. ( |3.43D has a nonzero macroscopic 
average, owing to the presence of a characteristic Pon- 
tryagin charge constructible from the Nj^, that feature 
in the energy shifts, and the magnetization direction. A 
second consequence of the SOI-generated carrier-energy 
shift ( 3.40| ) leads to the second contribution, cr^U. Due to 



the feedback of the (fast) carrier freedoms, which provide 
an effective p otenti al for the (slow) spin system, deter- 
mined by Eq. (3.40), the equilibrium probabilities of spin 
configurations having opposing Pancharatnam fluxes will 
no longer be equal. (For this contribution, which is re- 
lated not to dG/dsj but to G itself, there is no need 
to account for SOI-induced carrier-energy shifts in the 
current now being averaged over a nonsymmetric spin- 
configuration distribution.) A contribution with this ori- 
gin has also been considered in Ref. a 
arc of the same order of magnitude. 



(1) 

AH 



and 



,(2) 
"^AH 



I. Structure of the conducting network and the AHE 
resistivity 

We now consider the question of how the physics of 
elementary triads of Mn ions relates to the macroscopic 
properties of manganites. For hopping conductivity, the 
pathways taken by the current depends sensitively on 
the details of the configuration of the core spins, ow- 
ing to the sensitivity of the hopping amplitudes to the 
core-spin alignments. In particular, regions having local 
spin configurations that are aligned roughly opposite to 
the macroscopic magnetization of the sample tend to be 



avoided by the current. This fact renders rather subtle 
the procedure for averaging over equilibrium spin con- 
figurations, which must account for effects such as local 
spin correlations and excitations of various types (i.e. spin 
wave and topological excitations). 

En route to computing the AH resistivity, let us try 
to identify through which triads the AH current tends to 
fiow. Carriers tend to pass through regions of lower re- 
sistance. However, currents through regions with aligned 
spins do not lead to an AHE because the relevant Pan- 
charatnam flux through such regions is small. 



cluster R 




' — ^ — * 




/u 


...1 

* 1 ** 


y.-.^i^ 


/ ^ 


■H 

/■ 


1 




1 



cluster L 




hedgehog 



(b) 



FIG. 16. (a) Two clusters of sites, denoted L and R, in a 
fragment of the sublattice of Mn sites. Within a cluster the 
core spin orientations are roughly the same; but spins in dis- 
tinct clusters have significantly different orientations. Note 
that the boundary between such clusters is liable to contain 
spin configurations that resemble the single hedgehog config- 
uration shown in (b). The heavy dotted lines in (a) indicate 
a triad of spins that contribute to a hedgehog configuration. 
Sites in the upper part of the cell shown in (b) contribute to 
the conducting network and, correspondingly, the magneti- 
zation (per site) in this upper part of the cell (shown as an 
open-headed arrow) is roughly that of the sample (see the 
discussion in this section). 

Furthermore, in the ferromagnet-to-paramagnet transi- 
tion region, where the magnetization is only a small frac- 
tion of its saturation value, even those spins within the 
network responsible for longitudinal conductivity have 
orientations that are typically splayed, relative to one 
another. Nevertheless, if one were to consider only the 
spins in this network, their average magnetization would 
be larger than that of the entire sample, and the typi- 
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cal solid angles formed by triads of such spins would be 
rather small. Thus, any AHE originating in such tri- 
ads would be not the dominant contribution. Moreover, 
close to the metal-insulator transition, hopping paths 
through the sample that encounter spins having rather 
common orientations do not exist. In fact, there is ex- 
perimental evidence for the existence of moderately sized 
clusters of aligned spins coming from neutron scattering 
data at temperatures near the ferromagnet-to- 

paramagnet transition. These data indicates that clus- 
ters of spins do indeed exist in which spins are aligned 
over a length scale of roughly 10 A. Thus, it is reason- 
able to envision the magnetic configurations as compris- 
ing rather well oriented clusters of, say, twenty to thirty 
spins, with adjacent clusters having rather different spin 
orientations. Furthermore, theoretical estimates of the 
scale of magnetic fields relevant for colossal magnetore- 
sistance are consistent with the existence of such clusters. 
(Such clusters can be regarded as being large magnetic 
polarons |^^.) As the characteristic Zeeman temperature 
associated with a single spin-3/2 in a magnetic field of 
7T is 20 K, whereas the characteristic temperature asso- 
ciated with colossal magnetoresistance is roughly 200 K, 
one is lead to the view that clusters of correlated spins 
involve on the order of ten spins. 

Now consider two adjacent clusters of roughly aligned 
spins (e.g. clusters L and R in Fig. Even the con- 

ducting paths connecting these clusters contain bonds be- 
tween ions having significantly misaligned spins. These 
spins belong to regions of magnetic inhomogeneity, e.g., 
inhabited by hedgehog excitations, (an example of "lat- 
tice hedgehog" is shown in Fig. |l6|b) which defines the 
border of aligned clusters as shown in Fig. |l6|a. As men- 
tioned in the Introduction, hedgehog excitations are long- 
lived topological spin excitations, the existence of which 
is known to be important for ferromagnet-to-paramagnet 
transition ||5^ , |56|] , even in the three dimensional case. 

Within these regions of magnetic inhomogeneity there 
are triads of splayed spins. Let us now address the 
question: What is the characteristic splay? To answer 
this question, imagine dividing up the spins into those 
within clusters and those within the border regions. Even 
though the magnetization per spin in a typical cluster 
is greater than the sample-average magnetization per 
spin, the clusters are misoriented relative to one another. 
Thus, the contribution to the magnetization per spin of 
the sample coming from the spins in clusters is not guar- 
anteed to exceed the sample average and, indeed, it seems 
reasonable to assume that it is, in fact, not so different 
from the sample average. If so, then the magnetization 
per spin of the spins in the border region would also be 
roughly the sample average. We shall make the hypoth- 
esis that this is indeed the case. Then the magnetization 
of typical triads of ions in the border region can also 
roughly be taken to be the average sample magnetiza- 
tion. We shall denote by /3 the characteristic angle that 
spins in a triad form with the direction of the overall 
magnetization of the sample M: cos/3 ~ |M|/Ms, where 



Ms is the saturation magnetization of the sample. 

Let us now imagine how charge carriers move between 
clusters having quite different ion spin orientations. Such 
motion is necessary for the existence charge-carrier prop- 
agation between electrical contacts. We have sketched 
a typical instantaneous configuration of the spins in 
Fig. As one can see, the "upper" part of a inho- 
mogeneous region formed on a cubic lattice of sites can 
serve as a path for hopping between clusters. For reasons 
that we will now give, triads of sites within this (and sim- 
ilar) regions are effective contributors to the AHE: 

(i) The three spins in the triad have positive components 
along the direction of magnetization. (Recall that clus- 
ters with magnetization pointing against the majority 
tend to be avoided by the current.) This allows partic- 
ipation by these triad sites in the conducting network. 
If all three sites participate in the conducting network, 
the triad can be an effective source of an electromagnetic 
force that leads to a Hall effect. As we have discussed, 
the net magnetization of the triad is roughly that of the 
bulk, which makes it magnetically compatible with its 
neighbors. 

(ii) Typical triads of spins, being located as they are 
between several misoriented clusters are significantly 
splayed. Therefore, the solid angle formed by their spins 
(i.e. the area of the geodesic triangle formed by their ori- 
entations on the unit sphere, also known as the Pontrya- 
gin charge of the spin configuration) is substantial and, in 
fact, close to the maximum possible value given the con- 
straint that the triad magnetization (per spin) be compa- 
rable to the sample magnetization (per spin). Thus, we 
adopt as a caricature of the spin configuration in regions 
contributing to the Hall effect a picture of splayed triads 
of spins of known magnetization density, residing within 
tetrads of spins on a lattice plaquette, such as those de- 
picted by the lattice hedgehog configuration shown in 
Fig. This scheme, in which we consider tetrads of a 
given magnetization and then select triads of sites in a 
tetrad, seems to us appropriate, given the cubic struc- 
ture of the sublattice of Mn spins. However, alternative 
schemes (e.g. in which one considers triads themselves 
or other assemblies of splayed spins of a given magne- 
tization rather than tetrads, and chooses triads out of 
these assemblies) lead to almost identical results (e.g. for 
the scaling of the Hall resistivity, which we discuss in 
the present and following sections). This insensitivity to 
details is all the more natural, given that we are deal- 
ing with an atomically disordered system. We note that, 
because hedgehogs are topologically stable, they provide 
a mechanism by which the spin configuration can sus- 
tain strongly splayed regions that persist for durations 
much longer than the characteristic time for charge mo- 
tion. Hence, in their presence, on can accurately treat 
the charge motion as taking place with a background of 
inhomogeneous but essentially static spins, which renders 
consistent the adiabatic treatment of the dynamics of the 
spins. 



23 



Thus, we arrive at the notion of an optimal triad. An 
optimal triad is a triad of spins residing in a tetrad of 
four spins around a plaquette of the cubic sublattice of 
Mn ions and having the following properties: (i) The 
tetrad has the magnetization density of the bulk; and 
(ii) subject to this constraint, the spins of the tetrad 
are maximally splayed (i.e. subtend the maximal solid 
angle and, in fact, are configured symmetrically around 
a cone). Note that if the lattice were triangular then 
we would simply have adopted a definition of optimal- 
ity in terms of maximally splayed triads (rather than 
tetrads) of spins. As mentioned above, in disordered sys- 
tems (such as manganites), the distinctions engendered 
by such options are unlikely to have a strong impact on 
the physical consequences of the picture. 

The motion of charge carriers through optimal triads 
gives rise to the AHE. We note that these optimal tri- 
ads have properties quite different from those of optimal 
triads contributing to the OHE in doped, non-magnetic 
semiconductors |Q: in the OHE setting, only two sites 
in an optimal OHE triad are connected to the conduct- 
ing network, whereas in the optimal AHE triads all three 
triad sites participate in the network. (Indeed, if alter- 
natively, one of the sites is not a part of the conducting 
network, its spin must be roughly opposite that of the 
spins on the other two sites; such a configuration would 
yield only a small Pancharatnam phase.) 

The question may arise why triads within tetrads (and 
not, for instance, tetrads themselves) are considered to 
be the dominant source of the AHE. By a contribution 
from a tetrad we mean one involving four overlap in- 
tegrals. As these overlap integrals are small, owing to 
the localized character of the carrier wavefunctions and, 
thus, the contribution from tetrads is suppressed, rela- 
tive to that from triads. We note that distortions due to 
doping, particularly deviations of Mn-O-Mn bond angles 
from 180 degrees, facilitates tunneling between Mn ions 
via plaquette diagonals (see Fig. |l^). As was estimated 
in [Q, the amplitude of transfer along diagonals is 0.5 
of that between nearest neighbor Mn ions. Recent tight 
binding model parameterization of local density approx- 
imation (LDA) studies |80| show that hopping via diag- 
onals is even more important, and its amplitude is 0.82 
of transfer amplitude between nearest neighbor Mn ions. 

Having discussed the structure of resistive network, let 
us now calculate the longitudinal and Hall resistivities 
of manganites in the regime in which conductivity pro- 
ceeds by hopping (i.e. at temperatures above, as well as 
somewhat below, the ferromagnet-to-paramagnet tran- 
sition). The longitudinal hopping conductivity arising 
from phonon-assisted hops between sites i and j is given 

by 



{ne'd'/kBT)Wo^ cos\6/2), 



(3.44) 



[c.f. Eqs. ( pj] ) and ( |3.2b )], where d is the distance 
between sites. Here, Wq is the rate of phonon- 
assisted direct hops, and we have explicitly separated 



out the Anderson-Hasegawa factor cos^(0/2). Corre- 
spondingly, the (anomalous) Hall conductivity is given 
by [c.f. Eq.(pl)] 



{ne'd'/kBT)Wl' 



(3.45) 



where Wl^ is the rate of hopping between the two sites, 
and accounts for interference associated with both direct 
hopping and hopping via an intermediate state on a third 
site. Note that the quantity Wi includes three Andcrson- 
Hascga wa fa ctors [and so does the Hall conductivity given 
by Eq. ( 3.37| )]. The task of computing the Hall resistivity 
Pxy, which in the limit of Gxx ^ c^xy under consideration 
has the form 



Pxy 



'xy 



(3.46) 



then reduces to a determination of a ratio involving the 
direct and indirect hopping rates W^-' and Wl'' (as a 
function of the magnetization texture). As discussed in 
Sees. [HD and IHF, Wl'' involves two-phonon processes, 
whereas Wq involves only singie-phonon processes. Be- 
cause of this, dependence on electron-phonon coupling 
constant, phonon occupation numbers, and charge car- 
rier occupation numbers, cancels from the relevant ratio, 
W{' / {Wq , SO that this ratio can be written as 



W,/W^ = ahC/kBT, 



(3.47) 



where a is a numerical factor describing the multiplic- 
ity of the various carrier-phonon interference processes 
(see REf. [^ and Sec. [H D| ), the number of intermedi- 
ate sites, and the difference between nearest- and next- 
nearest-neighbor hopping amplitudes. We shall refer to 
the parameter which characterizes the difference be- 
tween the forward W{' (backward W(^) transition rates, 
as an asymmetry parameter. For the OHE, this asym- 
metry parameter is given by 



C oc sin(B • Q/(j)o) 



(3.48) 



where Q is the vector area defined in Eq. ( |3.42| ), as fol- 
lows from Eq. (3.3). 

For the AHE, it follows from Eqs.( |Oo| ) and ( |3.43| ) 
that the assymmetry parameter is given by 



C ^ 3[gjk ■ (rij X nfc)][ni • (na x n3)]/4, 



(3.49) 



where gjk are characteristic vectors arising in the hop- 
ping amplitude owing to the spin-orbit quantal phase; 
iij are unit vectors of the core spins in the triad, and 
ni • (n2 X na) is the volume of a parallelepiped defined 
by core-spin vectors, i.e., the Pontryagin charge qp. The 
anomalous Hall resistivity can be written in the simple 
form 



Pxy 



'^xyl O'xx 



1 f ahC 1 



ed2 cos" (0/2) 



(3.50) 
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The evaluation of Eq. (3.5C) reduces to a determination 
of 9, along with the products (rij x n^) and ni • (n2 x n^), 
which survive averaging over all possible triads. The 
dominant contribution to the average of these products 
arises from optimal spin configurations. 




(b) 



M 



FIG. 17. (a) Triad of Mn sites (1, 2 and 3) in the con- 
ducting network. Charge-carrier motion around triads such 
as this lead to the AHE. We compute the longitudinal and 
AH conductivities by relating them to the configuration of 
the spins (which we suppose to be optimal) and, hence, to 
the magnetization. 



Therefore, in line with properties (i) and (ii) of these 
configurations, consider a square lattice formed by the 
Mn ions in a plane perpendicular to m (as in, e.g., the 
top surface of the cube in Fig. p^). To ascertain the 
geometry of optimal triads, consider the spins at four 
sites of a plaquette belonging to the conducting network. 
Being optimally configured, these spins lie at equal sep- 
arations around a cone whose vertical angle 2/3 is given 
by 2cos"i[Af(iJ,T)/Msat]. (The angle between the alti- 
tude of the cone and any generator in the conical surace 
is p.) Let us now use this information to fix the vari- 
ous geometrical quantities that determine the longitudi- 
nal and Hall conductivities, the former being associated 
with pairs of sites and the latter with triads. We assert 
that to compute the contributions to these conductivi- 
ties we may consider sites in an optimal configuration, 
as these are characteristic of that regions of the sample 
that dominantly contribute. For the Hall conductivity, 
the reasons for accepting this assertion were discussed in 
the present subsection. As for the longitudinal resistiv- 
ity, the assertion is valid because this quantity is dom- 
inated by the most resistive regions of the conducting 
network, and these are expected to arise at the interface 
between clusters of aligned spins, i.e., in regions that are 
hedgehog-like. 

With this picture in mind, we now compute charac- 
teristic values for the geometrical quantities that feature 
in the longitudinal and Hall conductivities. Thus, we 
need the angle 9 between adjacent spins, the Pontryagin 
charge gp, and the products (rij x n/j), each of which is 
related to the vertical angle 2(3 by elementary geometry: 



2cos2(6l/2) 1 + cos2/3, 

qp = 2cos/3 sin^ /3, 
m • (rij X rifc) = sin^ /3. 



(3.51a) 
(3.51b) 
(3.51c) 



We now have to account for the fact that, in the 
hopping regime, the magnitude of the longitudinal (and 
anomalous Hall) resistivities depend on the probability 
that pairs and triads of ions are connected to the con- 
ducting network. We introduce a percolation factor P 
describing the connectivity of the pair to the conduct- 
ing network; for the AH conductivity the correspond- 
ing factor would be because both pairs of ions in a 
triad must, as discussed above, belong to the conduct- 
ing network. It is remarkable that, throughout the lo- 
calization regime, pxy is nevertheless determined by cur- 
rents formed in individual pairs and triads, because the 
factor s of P ca ncel in the expression for pxy given by 
Eqs. (|3l6|, |3^). Therefore, in so far as gp and the an- 
gles between neighboring spins can be directly related to 
m = M/Msnt = cos/3, the Hall resistivity p^y depends 
on H and T only through m{H, T), and is given by. 



Pxy 



m(l 

Pxy 



m 



2\2 



(1-Hm2)2 



(3.52) 



To determine the magnitude of the AHE, we first need 
to estimate the characteristic values of |gjfe| ~ g aris- 
ing fro m the spin-orbit interaction (SOI). As discussed 
in Sec. 



[HG 



the SOI term leads to a Dzyaloshinski- 
Moriya contribution to the eigenenergy of the carriers. A 
standard estimate |81| gives the characteristic values of 
|g| ~ 5 ~ Ze^ / 4:171 eC^o, where do is the radius of an Mn 
core d-state. While renormalization of carrier parameters 
in crystals may tend to increase |gjfe|, crystalline symme- 
try requires admixtures of core orbitals, which, in turn, 
are mixed with oxygen p-orbitals, with outer-shell wave- 
functions in order to have \gjk\ ^ . Such admixture is 
effectively generated by the non-collinearity of the Mn- 
0-Mn bonds that allows carrier hopping around triads 
(including jumps along plaquette diagonals). An estimate 
based on free electron parameters gives t; 5 x 10~^. 
(We note that band-structure calculation of the spin- 
orbit coupling constants is outside the scope of this pa- 
per). The characteristic strength of the Dzyaloshinski- 
Moriya terms is ~ gt^ ~ 0.02 meV, and the characteristic 
strength of the spin-orbit interaction is ^ eto ^0.1 meV, 
where e is the characteristic carrier energy. Not only 
these strengthss are much smaller than the characteristic 
double exchange energy, but they are even smaller than 
the magnitude of the direct antiferromagnetic Heisenberg 
exchange term. However, for the anomalous Hall effect 
in the localized regime, the Dzyaloshinski-Moriya terms 
are crucial, as we discussed in Sec. [II F. 



We now estimate the macroscopic longitudinal and 
Hall resistivities in the regime in which the conduct- 
ing network is fully connected, i.e., in the regime IV of 
Fig. |l^. By taking n = 5.6 x 102i cm^^ Wq ~ 2.5 x lO^^ 
s~^, and, from the magnetization data at T = 275 K (Fig. 



25 



1), COS P = 0.6, we obtain p^x — 1 cm which coincides 
with the value of the experimentally observed resistivity 
for LPMO (see Fig. |). The AHE contribution to the 
Hall resistivity, assuming a numerical factor a of 2.5, is 
then pxy ~ —0.5 /ifi cm, in agreement with the experi- 
mentally observed LPMO Hall resistivity at the same T 
(Fig. 1^). The equivalent expression for the hopping Hall 
resistance in the Holstein mechanism is defined by the 
asymmetry parameter C — cos^(0/2) cos/3sin(B • Q/i/io) 
and, at -B = 1 T, is an order of magnitude smaller than 
the AHE. We expect the macroscopic hopping AH and 
OH effects to have the same sign, opposite to that of the 
OHE in the metallic regime. 



In the next section. Sec. IV, we shall compare the re- 



sults for the Hall resistivity with the experimental data. 
As we shall see, the picture for the core spin configura- 
tions developed above, which include clusters of oriented 
spins and hedgehog-configured spins, allows us to explain 
not only the AHE, but ferromagnet-to-paramagnet and 
metal-insulator transitions in manganites, and provide 
a quantitative explanation of the magnitude of charac- 
teristic magnetic field that result in colossal magnetore- 
sistance. The notion of an optimal triad enables us to 
fit the experimental data for the AHE to a functional 
depe nden ce of the resistivity on magnetization given by 
Eq. (3.52). The agreement of the hopping picture and 
experimental data in the transitional region is remark- 
able. 

As we have mentioned above, the structure of the con- 
ducting network leading to the ordinary Hall effect in 
disordered doped semiconductors and the averaging pro- 
cedures in these systems are still controversial |7^j7J . In 
contrast to disordered doped semiconductors, mangan- 
ites turn out to be systems in which the ability to tune 
average magnetization allows one to tune optimal triads, 
whose solid angles (the Pancharatnam phases) determine 
the AHE. The magnetization in manganites, therefore, 
serves as a scaling variable that has no analog in OHE in 
nonmagnetic disordered systems, and provides a check 
on our understanding of the conduction network. We 
note that the presence or absence of small polarons in 
the system does not change the scaling of the AHE resis- 
tivity, because, as follow s from studies of polaronic trans- 
port in [imiJa^,!?!! , Eq. |3.47| also holds for small polaron 
hopping. 



are swept from the sample before saturating and follow- 
ing the general trend. At the lowest temperatures, the 
metallic OHE appears as a positive contribution at con- 
stant magnetization. As for the magnitude of the Hall 
resistivity, 
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FIG. 18. Hall resistivity p^y of LCMO versus reduced mag- 
netization M/Msax for the data shown in Fig. 4. Note the scal- 
ing behavior, i.e., the extent to which p^y can be regarded as 
depending on T and B solely through M/Msat. The solid line 
is a fit to Eq. (3.52) with p^y = -6.2 pQ,cra. 



for LPMO the solid curve in Fig. |T| follows Eq. ( |3.52| ) 
with p^y — —4.7 pVl cm is consistent with the estimates 
of pxx and pxy given above. 



a -0.4 

of -0.6 
-0.8 
1.0 



lao K.aiD K 




LPMO 



0.0 0.2 



0.4 0.6 0.8 

/M . 



1.0 



IV. HALL RESISTIVITY: COMPARISON OF 
THEORY AND EXPERIMENT 



The scaling of the Hall resistivity is shown in Figs. |T^, 
19| and in which the data shown in Figs. ^, |5| and 
Irrespectively, are replotted as a function of M/Ms^t- 
At and above Tc the data fall on a smooth curve that 
reaches an extremum at M/Mgat — 0.4 for LSMO and 
LPMO and at M/M^at ^ 0.35 for LCMO. Below Tc the 
data first change rapidly with magnetization as domains 



FIG. 19. Hall resistivity p^y of LPMO versus reduced mag- 
netization M/Msat for the data shown in Fig. 5. Note the scal- 
ing behavior, i.e., the extent to which pxy can be regarded as 
depending on T and B solely through M/Msat. The solid line 
is a fit to Eq. (3.52) with p°y = -4.7 pQ.cm. 

Down to 285 K, which is the Curie temperature deter- 
mined by scahng analysis, Eq. ( 3.50| ) describes the data 
for LPMO reasonably well. In addition, the extremum 
found from this equation is located at M/Msat = cos (3 « 
0.35, close to the experimental extremum. 



26 



In LCMO and LSMO, the agreement between theoret- 
ical and experimental results is good as well. We note 
that below Tc, the longitudinal resistivity is metallic and 
no longer dominated by magnetic disorder. We have not 
expected an agreement between theory and experiment 
in this range of temperatures, but in LCMO and LSMO 
scaling persists at temperatures below Tc, with notable 
exception of the range of low temperatures and magneti- 
zations close to saturation value, where the ordinary Hall 
effect manifests itself. We note in this regard that below 
Tc, local spin arrangements still can still dominate the 
AHE via asymmetric scattering or side jumps. 



explains the results quite well. The anomalous Hall ef- 
fect arises due to interference between direct hopping be- 
tween two sites and hopping via a third site, with the 
quantal phase provided by topologically nontrivial con- 
figurations of Mn ion core spins in the presence of strong 
Hund's rule coupling. These force the hopping charge 
carrier to follow the local spin texture, with the average 
quantal phase arising due to local Pancharatnam phases 
and Dzyaloshinski-Moriya spin-orbit interactions. Below 
the transition temperature, the AHE competes with the 
OHE as long-range magnetic order and, presumably, an 
infinite percolating metallic cluster, develops. 
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FIG. 20. Hall resistivity p^y of LSMO versus reduced mag- 
netization M/Msax for the data shown in Fig. 6. Note the 
scaling behavior, i.e., the extent to which p^y can be regarded 
as depending on T and B solely through M/M^ax- The solid 
line is a fit to Eq. (3.52) with p'^y = -3.4 pTlcm.. 



The numerator of Eq.( 3.5^ ), m(l — m^)^, which is char- 
acteristic for the behavior of Uxy alone, has an extremum 
at m = 1/a/5 ~ 0.45 as shown by the dashed line in 

Fig. m. 

The broader maximum in the data suggest a shift to- 
ward a hopping model for p^x and pxy as the sample is 
warmed through the metal-insulator transition. 



V. CONCLUSIONS 



Our investigation of the Hall resistivity, the longitudi- 
nal resistivity, and the magnetization in single crystals of 
three different manganite compounds suggests that near 
and somewhat above the ferromagnet-to-paramagnet 
transition temperature, transport properties are deter- 
mined by charge carrier hopping between localized states. 
We find both theoretically and experimentally that the 
Hall resistivity is solely determined by the sample mag- 
netization (M) near and somewhat above the transition 
temperature. A microscopic model for the anomalous 
Hall effect based on the Holstein picture of the ordinary 
Hall effect in the hopping regime has been proposed and 
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